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ABS eve 


Relative vorticity fields calculated from the U. S. Navy operational Global Band 
Analysis are used to relate synoptic and storm parameters to the track of tropical 
cyclones in the western North Pacific Ocean. In this preliminary study, synoptic 
patterns are developed, described and discussed from the perspective of a pattern 
recognition technique to assist the forecasters at the Joint Typhoon Warning Center, 
Guam. The focus is on track turning motions to the left and right of the persistence 
track and on trving to accurately predict the point of the turn or recurvature in relation 
to the time evolution of the vorticity patterns. [he developmental sf; ¢ Of Segmmme 
indicates that there is potential for using synoptic patterns in the Global Band Analysis 


to guide the selection of the appropriate track aid in the 48-60 hour time range. 
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A. BACKGROUND 

The successful forecasting of tropical cyclone motion is one of the most difficult 
problems in tropical meteorology. The training process for the operational forecaster 
includes a process of integrating knowledge contained in the literature and learning 
from experienced forecasters. One aspect of the training is the development of an 
ability to recognize significant patterns in the synoptic analyses, satellite imagerv and 
forecast fields. The better forecasters are able to recognize the synoptic patterns that 
indicate that significant changes will occur (or not occur) in the track of the storm 
during the forecast period, or are able to recognize the synoptic patterns earlicr and 
more often than the less successful forecasters. 

A significant problem at the Joint Typhoon Warning Center (JTWC), Guam, is 
the short tour length for Air Force and Naval officers assigned as forecasters. Nearly 
half of the forecaster’s tour at JTWC is spent learning the job. The actual time 
required to complete the basic training phase depends on the abilitv and experience of 
the forecaster trainee and the number and type of storms that occur during the training 
period (McLawhorn, 1984). The short tour length at JTWC 1s not the only obstacle to 
fieemne the COMSEVENTHPFLT requirement to reduce forecast errors to 50, 100 and 
150 n mi for the 24, 48 and 72 h forecast periods respectively. The National Hurricane 
Center (NHC), Miami experienees forecast errors of magnitudes similar to those at 
JTWC without having to deal with the problem of forecasters subject to 15 to 24 
month rotation periods as in the military. JIWC and NHC each have several 
forecasting aids available for the development of the operational forecasts. Both centers 
are Continually testing new methods developed by the research community, as well as 
making improvements to current operational models. 

It 1s sometimes stated that the forecaster does not need more forecast aids. 
Rather, the need is to devclop the ability to better utilize the schemes that are available 
in the forecast centers (Neumann and Pelissicr, 1981; Tsui, 1984; Elsberrv and Peak, 
1980). Each of the forecast aids has strengths and weaknesses that depend on a 
complex combination of the storm size, speed, intensity, location and the synoptic 


features in the basin. The features of importance in the western North Pacific Ocean 


region include the subtropical ridge, long-wave troughs. short-wave troughs. the 
Northeast and Southwest monsoons, the upper-level flow patterns and features, and 
Other tropical svstems. Lach of these featiires is Gispiiec 900 SOiieire .t0 ie 
analysis and prognostic fields tha are used by the forecaster, who must decide whieiiigm 
the forecast aids is the best at depicting the interaction of the environmental [remiss 
the storm. In many cases, the more subtle features that are undetected by the 
forecaster under pressure to issue the forecast are later “discovered” in the post-storm 
analvsis. Successful forecasts both incorporate a maximum of information into a 
warning track and help the operational customer to make prudent and timely decisions 
to save life and property at a minimal monetary cost. 

Consider the set of forecast aids available for the [2 UTC 19 August 1983 
warning on Tropical Storm Dom (Fig. 1.1). This ts not an atypical pattern and makes 
the selection of an accurate forecast track very difficult. What is required is a method 
that will incorporate the storm and synoptic parameters into a pattern or set of 
patterns that are easily recognizable under the pressures of getting the forecast out. 
This will assist in the selection of the primary forecast aid and simplify the track 
forecasting dilemma for the operational forecaster. This feasibility study will seek the 


synoptic pattern relationships involved in the recurvature or nonrecurvature forecast. 


B. RECENT STUDIES IN SUPPORT OF THE JOINT TYPHOON WARNING 
CENTER 


One attempt at developing a technique that assists in the interpretation of the 
environmental fields is the CYCLOPS Objective Steering Model Output Statistics 
(COSMOS) forecast aid (Allen, 1984), which uses the output from the CYCLOne 
Prediction System (CYCLOPS). The success of this technique in predictimemanie 
unusual track of Tropical Storm Dom is shown in Fig. 1.1. COSMOS is a model 
output statistics technique applied to the CYCLOPS geostrophic steering flow aid, 
which is the most recent version of the HATRACK/MOHATT steering flow model 
(Renard, 1968; Renard and Levings, 1969; Renard er al., 1970; Renard ef Gigi 
Renard et al., 1973). The COSMOS program compares the current flow pattern to a 
set of seven typical patterns and calculates a forecast track that incorporates the 850, 
700 and 500 mb steering flow outputs from the CYCLOPS model. The output of a 
single track (versus separate tracks at each level) simplifies the forecaster’s decision 
problem, incorporates past pattern information and allows the displav of the 


CYCLOPS components that are used in the technique. 
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Another recent attempt to incorporate the environmental circulation pattern into 
a forecast aid is the Coiorado State University Model (CSUM) developed bv 
\latsumoto (1984). This staustical technique categorizes the track forecast problem 
Into one of three classes depending on the storm center location relative to thie 
sudtropical ridge at 300 mb. The three classes are defined as under the ridge. on the 
ridge and above the ridge. The predicted track is divided into three 24-hour periods 
and uses the prognostic fields as if thev are analvses at the valid ume. 

A series of investigations at the Naval Postgraduate School over the past three 
Vears has attempted to improve the use of forecast guidance at JTWC. Currv (1985) 


developed an objective technique to improve the initial location of the storm center 


that is used for the forecast aids. Wilson (1984), Schott (1985), Williams (1986), Jones 
(1986), Meanor (1987) and Wenrger (1987) have studied various combinations of storm 
and SViODtic puranv@ters to Cee cian information for possible inclusion in w@ decisionstiee 
scheme for selecting the optimum track forecast guidanec fereaspartieular loreemme 
scenario (Peak and Elsberrv, 1987). These studtes have inereased the knowledgesbage 
for the forecaster and researcher, but have not produced a straight-forward pattern that 
is easily recognizeable by the forecaster. 

Sandgathe (1987) has categorized and described the various storm-synoptic 
relationships that appear to have the greatest influence on the tropical cyclone tracks 
in the western North Pacific Ocean. He has also reviewed the 1982-1985 storm seasons 
and listed the interactions that affected the tracks of selected storms. A detailed 
description of the interactions and the forecast problems that arise due to the 
occurrence of the interaction, or a multiple interaction, are given in the report. A 
description of the interactions studied in this investigation are given in the section on 


the data set development in Chapter IT. 


C. OBJECTIVES 

This study is a preliminary investigation of the storm and synoptic parameters 
related to tropical cyclone recurvature. The goal is to develop a way to display and 
analyze the parameters in a manner that 1s useful to the operational forecaster. The 
recurvature problem is chosen because of the difficulty experienced in accurately 
forecasting the recurvature location and time. An inconsistent performance of this 
critical forecast task often results in 72 h position errors in excess of 1000 n mi. These 
errors far exceed the COMSEVENTHELT requirements and create major problems for 
operational commanders tasked with making decisions to implement typhoon evasion 
and preparation plans. Reduction of these track errors to levels comparable to the 
current performance for straight-moving storms would provide a credible and useful 
product to the operational decision makers. 

Several approaches to the problem have been used to discover a more optimal 
method of identifying and displaying the storm and synoptic parameters that determine 
the track of the tropical cyclone. The data fields, data retrieval and display are 
described in Chapter II, analysis techniques and results are presented in Chapter ITI. 
Chapter IV discusses the conclusions and possible avenues for further development and 


study. 
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ti. DATA BASE 


A. GLOBAL BAND ANALYSIS AND STORM DATA 

The wind fields in this study are taken from the United States Navy Fleet - 
Numerical Oceanography Center (FNOC) Global Band Analvsis (GBA). The GBA 1s 
the analysis used to initialize the Navy’s two dynamic tropical cyclone track forecast 
models, the Nested Tropical Cyclone Model (NICM) and the One-way Influence 
Tropical Cyclone Model (OTCM). The GBA are produced at 00 and 12 UTC each day 
on a 49 x 144 grid with grid spacing of 2.5 degrees of longitude on a Mercator 
Peoecton., le erid extends from: 40.956°S to 39.745°N and has 360 deerces 
longitudinal coverage. The grid spacing decreases toward the poles as the meridians 
Converge. At the standard latitude of the grid (22.5°N or S$), the grid distance is 257 
km. | 

A complete description of the GBA can be found in the Numerical 
Environmental Products Manual (United States Naval Weather Service, 1975) and a 
review of recent developments in the Navy's tropical analysis is given in Elsberry and 
Fiorino (1985). The analysis is done at the surface, 700 mb, 400 mb, 250 mb and 200 
mb, Thermal analyses at intermediate levels provide data for the vertical coupling 
using the thermal wind relationship. The operational centers receive intermediate level 
fields calculated from analyzed levels, but these fields are not archived. When a 
tropical cyclone is present, the storm is depicted by eight wind vectors bogussed into 
the surface field around the position of the center at chart time. Data used in the GBA 
are from the surface observation network, rawindsondes, pibals, aireps and cloud 
motion vectors. If no observations are available for the current analvsis, a blend of the 
12 hour old analysis with a 5% reversion to climatology is used. 

Data for tropical cyclones in the western North Pacific basin are available on an 
annual basis from the JTWC’s Annual Tropical Cyclone Report (ATCR). The report 
provides a narrative description of each storm with the storm best track, satellite 
pictures, significant synoptic patterns and occasionally storm damage photographs. 
The data annex of the report gives six-hourly best track, warning and forecast 
positions; forecast error statistics; satellite, aircraft and radar fixes; best track intensity, 
aircraft wind reports and satellite intensity estimates. In this study, ATCR’s for 1982, 


1983 and 1984 were used for case selection and analysis. 


B. STORM SELECTION 

Super Typhoon Marge (1SW) 1983 was chosen as the first stormy to be Sttigiem 
due to the difficulty experienced bv JTWC in accurately forecastilig the Storie mee 
throughout the life of the system (Fig. 2.1}. Marge was categorized bY Santis 
(1987. p.19) as a step-recurver that later underwent extratropical transition. A step- 
recurver is a system that is influenced by the passage of mid-latitude troughs to the 
north of the storm. The step phase of the track takes place as the storm slows and 
moves northward under the influence of a passing trough, but does not recurve. As the 
trough passes or weakens, the subtropical ridge builds westward and the storm turns to 
the left and tracks to the northwest or west under the ridge. Marge stepped twice 
‘before the eventual recurvature. Each time the step started, large forecast errors Were 
recorded as forecasters predicted recurvature. When the storm finally did recurve, the 
forecasters unfortunately continued a track toward the northwest, which again added 


to the overall error. 


TABLE 1 


TRACK FORECAST ERRORS (N MI) IN 
SUPER TYPHOON MARGE COMPARED TO 1983 OVERALL ERGs 


Warning 24h 48h JQ) 
Marge BY ea 484 HS 
1983 16 a) oo 405 
Marge/ 1983 1.18 1.63 1.87 1.86 





Inconsistent forecast tracks as found in the Marge case can be responsible for 
premature evacuations of aircraft and ships and the early or unnecessary setting of 
tvphoon conditions. Both actions entail major expenditures in terms of both man-hours 
and money. Table 1 is a summary of the average track errors for Marge compared 
with the overall errors for 1983. Although the errors and the ratio to the 1983 average 
are large, there were even larger errors on other 1983 storms in all categories exceputhe 
72 Vi toneeast. 

A set of storms was then compiled to test the preliminary results based on 


Marge. The goal was to develop a set of storms that would allow the study of varying 
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Fig. 2.La Super Typhoon Marge best track (dots) and forecast tracks (stars) during 
the first step. Labels indicate date and hour (02 00 means 00 UTC 02 Nov). 
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Fig. 2.16 Super Typhoon Marge best track and forecast tracks 
during the second step. 





P20 


© 
0700 


¥ O8 00 
06 12 





0712 ya O TOO 





Fig. 2.lc Super Typhoon Marge best track and forecast tracks 
during the recurvature. 
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storm tracks to test the pattern relationships between the svnoptic ficld and the storm 


envelope and to include as anv storm-svnoptic interactions as possible. 


Sandeathe (1957) has descertved the traca [forecusang provleny 11 ie Ves 


North Pacific Occan in terms of the storm-synoptic interactions, bis list shia 


ranked by need for improved foreeaster understanding, includes: 


(111) 


(Vii) 


Multiple cyclone interaction -- the flow created or modified by one cyclone 
influences the track of the second. 

Cyclone - midlatitude trough interaction -- the trough influences the 
cyclone track and causes a movement to the north and then recurvature if 
the trough influence is strong enough, or a step pattern if the trough 
influence is overcome bv a strengthening of the subtropical ridge. Looping 
may also result from this interaction. 

Cyclone - subtropical ridge interaction -- the shape and extent of the ridge 
dictates the path of the cyclone to the west or in a gradual recurvature 
around the end of the ridge. The ridge also influences the step and loop 
patterns. 

Extratropical transition -- the warm core cvclone transforms into a hybrid 
or cold-core system. 

Terrain interaction -- the track and intensity of the cyclone are influenced 
by interaction of the steering flow or the storm circulation with the various 
major islands and coastal areas in the Pacific basin that have significant 
orographic features. 

Monsoon surge interaction -- the tracks of some tropical storms show 
diurnal variations due to the diurnal fluctuations in the monsoon flow 
during the Southwest monsoon, which sometimes causes the storms to 
become quasi-stationary in the South China Sea and Philippine Sea. The 
Northeast monsoon also hays a strong influence on the track of late season 
tropical storms, and may lead to erratic tracks and large forecast errors. 
TUTT or upper low interaction -- the influence of the Tropical Upper 


Tropospheric Trough (TUTT) lows on cyclone motion. 
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The 12 storms in the data set and the interaction categories according to Sandgathe are 
listed below (numerals following the storm name indicate the storni number and vear 


and the Roman numerals correspond to the interaction list above): 


i ©6=l Besse T1782) (11,111 

(2) Gordon (16W;/82)  (1,11,11) 

(3) Wayne (04W/83)  (v,vu) 

(4) Abby (05W/83) (1,11,111,1V,V) 
(5) Forrest (11W/83) = (a,1V,vu1) 


(6) Ida (14W/83) (11,1V, Vil) 
arm Lex (17We83) (11.V1) 
($)} Marge (18W,83)  {itiv) 
ee) tke (15. WS4) (111) 


(10) Thad (24W/84) (11, V1) 

(11) Vanessa (25W/84) (,u1,1v) 

(12) Clara (29W/84) (1,11) 
In addition to Super Typhoon Marge, the set includes five storms that have the step 
pattern, five that recurved without first stepping and two storms that had no apparent 


interaction with the mid-latitude troughs and had straight tracks. 


C. DATA RETRIEVAL 

Wind fields are extracted from the GBA archives at 00 and 12 UTC for each of 
the storms studied. The fields include the zonal (u) and meridional (v) wind 
components at the surface, 700 mb, 400 mb and 250 mb. 

The data tapes are processed in a series of steps to allow the display of the wind 
data and derived fields in the Interactive Digital Environmental Analvsis Laboratory 
(IDEA Lab) at the Naval Postgraduate School. The IDEA Lab consists of a 
VAX8200 mini-computer and two nucro-VAX computers linked to various display 
terminals and input/output devices. The system allows interactive display of satellite, 
observational and gridded data as well as a variety of data processing capabilities 
(Wash et al., 1987). The first step in the process is to match the fields and storm 
tracks. A list is compiled using the ATCR best track data for each storm and 1s 
compared to the fields available on the data tape. A full listing of the fields available 
can be produced using the Global Band File Identification (GBFID) program found in 
Appendix A. 


I] 


The next step is to extract the fields required for study fron them. 4) 
Environmental Display Network (NEDN) format tapesmandmereate” a iane in /voiee 
format for use in the IDEA Lab. This is accomplished using the program GETGBA 
listed in Appendix B. The data section of the program niust be cureiully prepaneuame 
avoid problems due to missing fields in the GBA archives. If a complete set of fields 
for a date and time are missing from the set, that date-time code entry must be omitted 
from the GETGBA data file. If only a partial set of levels or component fields is 
missing, the date-time code information should be entered, which allows the 
subroutines to extract the fields that are present. The number of records to be skipped 
to get to the desired fields must also be accurate (see documentation in Appendix B). 

The ASC IT tape is used to create the data base in the IDEA Lab. The creation 
of the files for display and manipulation of the data is a two part process. First. the 
tape data 1s read to storage on a disk using the program: Global Band Tape To Disk 
(GBTTOD) found in Appendix C. The conversion of the data to a format that 1s 
compatible with the graphics software is accomplished using the FILE CONVERTER 
program in Appendix D. This program changes the GBA grid to a 2.5 degree latitude 
by 2.5 degree longitude grid using an interpolation subroutine, and inverts the grid 
array to make it compatible with the GEneral Meteorological (GEM) software 
(developed at NASA’s Goddard Space Flight Center) used in the IDEA Lab. Oncesthe 
data are written to a GEMFILE, the fields can be displayed and manipulated using the 
programs GEMPAK and GEMPLT. 


D. GEMPAK AND GEMPLT PROGRAMS 
The GEMPAK and GEMPLT software are used to display the data in a variety 

of ways including tabular data, grid point plots of data, contour plots, vector plots, 
Streamlines and vertical profiles. The gridded data can also be manipulated with a 
variety of Grid Diagnostic Programs to display numerous derived fields including: 

a) Sums (ADD) 

be»  Dikiterence (SUB) 

c) Laplacian (LAP) 

d) Laver Average (LAV) 

Cc) Time Derivative (DDT ) 

f) Time Difference (TDF) 

g) Time Average (TAV) 

h)  Drtvergence (DIV) 

1) Relative Vorticity (VOR) 


12 - 


}) Absolute Vorticity (AVOR) 

kK) Gradients (GID) 
Pees iitber O01 Ofer lunctions found in the GENIPAK User's Guide (desJardin ev a. 
Pro ssomero! the di@enostic grids can be created “omthe flv" and be directlv displaved 
Or sent to a hard copy output device. More complicated functions require that 
intermediate grids be created before further processing. In either case, an option exists 
to save the diagnostic grids in the same GEM file as the original storm data in date- 
time order. 

The display terminals give multi-color depictions of the fields including 
geography of various resolutions and projections. Fields can be overlaved on the 
screen for analysis and can be reproduced in hard copv on laser printers or full color, 
Seeae ¥IERSATEC plotter. 
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til, DATA ANALYSIS 


A. ANALYSIS PARAMETERS 

Analysis of the fields is based on significant storm and synoptic parameters and 
the changes with time in the synoptic fields. The storm parameters used in the study 
are those used bv the operational forecaster to request and initialize the forecast aids 
generation programs at FNOC. The parameters used in the study (Table 2) differ 
slightly from the combined aids request (CARQ) as the storm position information 1s 
reduced from 48 to 36 h and the intensity and size data are increased to 36 h of history: 
The storm size information 1s taken from the warning messages using the wind radu for 
30, 50 and 100 kt. The values relating to the outer closed isobar as required in the 


forecast aids request were unavailable. 


TAB DE? 
STORM PARAMETERS 


. Month of year, dav of month 
. Position 
a. Latitude, longitude 


b. Fix platform (aircraft, satellite, radar, synoptic or extrapolation) 


c. Position and direction change - movement past 12, 24, 36h 


. Intensity 

a. Surface wind (kt), MSLP (mb) 

b. 700 mb winds (kt), height (m) 

c. Change in above parameters in the last 12, 24, 36 i: 
eS LOMIESIZE 

a. Radi of 100, 50 and 30 kt winds 


b. Change in radu in the last 12, 24, 50h 
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The synoptic parameters (Table 3) include those required for the CARQ as well 
aseicatunes seu il Other “forecast techniques or deemed significant through 
Perelomimentar research findings, Tne position and intensity of the subtropical ridve 
and nud-latitude trough are included as they are part of the CARQ values, and the 
history change values are extended to 36 hours. The polar and subtropical jets are 
considered in light of the Typhoon Acceleration Prediction Technique (TAPT) 
developed by Weir (1982). The positions of the short wave trough north of the 
tropical cyclone were added to the list after the initial analysis had begun, as small 
deviations in the track were often found to coincide with the short-wave trough 
passage. As in the storm parameters, the changes in the synoptic parameters are 


recorded for the past 36 h. 


B. STREANILINE FIELD ANALYSIS 
The initial analysis of the data involved a detailed study of Super Typhoon 

Marge using the streamline version of the GBA wind fields (as used at JTWC) and the 
National Meteorological Center (NMC) Tropical-Mercator surface and upper-air 
analyses. This approach proved difficult for several reasons. The NMC analyses were 
not available for all time periods, and the mid-tropospheric chart was for the 500 mb 
vice the 400 mb level in the GBA series. This led to problems in developing a 
continuous history and deternuning accurate height center locations corresponding to 
the wind field features. The GBA streamline output presented inconsistencies in the 
depiction of some small centers, and often showed cyclones and anticyclones with 
reverse flow patterns. While the bearings from the storm center to the synoptic 
features could be determined in a consistent manner, the range determinations were 
often suspect due to scale changes across the Mercator grid. After completing the 
analysis of Marge using this approach, several conclusions were made: 

a. [he process requires too much time to be useful in a forecast situation; 

b. The results are not clear due to the large number of variables involved; and 


c. Inconsistencies fron: chart to chart reduce the value of the information. 


C. VORTICITY FIELD ANALYSIS 
1. Vorticity Field Presentations 
The storm and synoptic parameters are important and require a presentation 
format that will incorporate the variations in the parameters into a more easilv 


interpreted form. This idea led to the use of the vorticity fields (Fig. 3.1), which are 
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SYNOPTIC PAR SIEGES | 


1. Subtropical ridge (Surface, 700 and 400 mb) 
a. range and bearing to the maximum height of the ridge 
b. range and bearing to the lowest height or break 
c. height magnitude range and extent of the ridge 
2. Long-wave trough (Surface, 700 and 400 mb) 
a. range and bearing to the trough axis at 35 deg N to the west 
b. height of the trough 
c. closing rate of trough to the storm center 
3. Short-wave trough 
a. range and bearing to trough to the NW-NE of the storm center 
b. closing or opening rate 
4. Upper-level trough 
a. TUTT location 
b. closed cell locations and intensity 
5. Subtropical and polar jets 
a. range and bearing to the jet max NW-NE of the storm center 
b. TAPT pattern type (Weir, 1982) 
6. Second, third or fourth storm in the basin 
a. bearing and range to each storm 
b. closing rates 
Cc. pressure and wind speed differentials 
7. Topographic influence 
a. range and bearing to the nearest influential land feature 
b. flow direction and speed over the feature 


c. closing or opening rate between the storm and feature 


presently not available to the JTWC forecaster. The bogussed storm and the analyzed 


features in the GBA wind field are displayed well in the relative and absolute vorticity 
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fields using the diagnostic capabilities of the GEMPAK software in the [DEA Lab. 
The conversion to vorticity climimated the reverse flow around small features that 


1 


aa ocean themsimeamlinessind allowed better analvsis of the details in the muaycr 
Poe comeliiemsiibtropical dee in the vorticity fields sives details of the relative 
strength of different areas of the ridge that are not represented well in the streamlme 
chart. The intensity of the synoptic troughs is also better shown in the contoured 
vorticity chart. The packing of the contours around the low centers shows the relative 
strengths of the centers as well as the changes in intensity when viewed in a time 
sequence. Thus, the vorticity presentation conveys important information to the 
forecaster at a glance without requiring lengthy data tabulation. The smaller scale 
storm and synoptic features appear to be integrated into the vorticity fields and 
Peeeemicd thirouel) the variations found in the contour patterns between the major ridge 
and trough depictions. 

It is hypothesized that the presence and location of a weakness in the ridge 1s 
very important when trying to predict the longitude of a turning or recurvature event. 
Perhaps an improvement in the prediction of the recurvature point could result from 
using the vorticity charts in conjunction with the objective track forecast aids. 

DeMaria (1985) suggests that the storm motion can be related to the absolute 
vorticity gradient of the steering current with motion components in the direction of 
the gradient and cross-gradient. Thus the analysis of the vorticity fields was done 
initially using the gradients of relative and absolute vorticity. The model results of 
DeMaria were not evident in the operational fields, perhaps because of errors and 
inconsistencies in the operationally-anatyzed fields that are then magnified bv taking 
gradients. Since the error characteristics of these vorticitv fields are not well known, 
the forecaster must determine whether a vorticity feature represents a real circulation 
or is caused by data inconsistencies. 

The relative and absolute vorticity fields contain several dominant patterns 
that are found to vary gradually in time as the interactions among the tropical storm 
and the adjacent synoptic features strengthen and weaken with time. The relationship 
of the cyclone track to the synoptic patterns revealed in this study (Fig. 3.2) are: 

(1) CHANNEL -- linear pattern with relatively straight relative vorticity 
contours on each side of the storm that results in a “cork in a streanv type 
pattern (Fig. 3.2a). That is, the past and future storm positions tend to lie 


Within the same “channel” of relative vorticity. 
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(11) 


(111) 


(1) 


NODE -- bulges in the contour(s) around the storm center that extend 
outward from the average circuinference of the contour around the center, 
or extend northward from the average latitude of a zonally-oriented 
contour to the north of the center (Fig. 3.26). Pmmary nodes are defined to 
be along the major axis of an elliptical vorticity contour that encloses the 
storm, or along the persistence track. That 1s, the tropical cyclone appears 
to track toward a primary node. However, secondary nodes are also 
menitored because they may become the major axis of the elliptical center, 
and thus be related to the subsequent storm direction changes with time. 
CLOSED CONTOURS -- vorticity (positive and negative) centers adjacent 
to the storm center that potentially mav influence the long-term track (Fig. 
3.26). 

PINCHED or BROKEN RIDGE -- areas of contour deformation in the 
ridge north of the center that indicate weaknesses in the ridge (Fig. 3.2d). 
STORM-STORM INTERACTION -- pattern transition from independent 
closed contours enclosing separate storm systems to a continuous Vorticity 
contour enclosing two or more storm centers (Fig. 3.2e). 

RESULTANT VECTOR PATTERN -- a closed vorticity contour vaitiva 
primarv and a secondarv node with axes approximately 80-120 deg apart 
(Fig. 3.2f). In some cases, a resultant vector from these two no@esse@aiere 
constructed that appears to have very good correlation with the track at 
60-72 h. 


2. Proposed Synoptic Analysis Technique 


The proposed analysis technique involves a recognition of the dominant 


svnoptic pattern from the above list. Changes in the synoptic pattern with time are 


then recorded and are related to expected changes in the storm track. The Omeniaimigm 


and deformations in the zero contour of relative vorticity is the kev to the analysis, as 


it is common to all plots and usually encloses the storm center. Based on George 
(1975), Chan and Grav (1952). Chan (1984) and Chan (1985), vorticity fields at 700 


mb, 400 mb and the laver average from 700-400 mb are examined to define synoptic 


pattern relationships to the track of the storm. The layer average is more consistent in 


time than the 700 and 400 mb fields alone. The laver average gives a smoother contour 


field that seems to change more conservatively through time than either of the 


individual levels, as suggested by Chan and Gray (1982). The relative vorticity fields 
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Fig. 3.2a Relative vorticity contour - channel! pattern. 
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Fig. 3.2c Closed relative vorticity contour pattern. 
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Fig. 3.2f Dual node-derived resultant vector pattern. 
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are chosen because the patterns of interaction among the storm and synoptic features 
are clearer with the earth's vorticitv removed. 

In the development phase of the research, a set of relative vorticity charts was 
plotted for the duration of Super Typhoon Marge. The best track for the 12 h old 
through the 72 h future positions were plotted on the charts to allow comparison of 
the vorticity pattern and the storm track. Unfortunately, the one chart set that was 
not available (12 UTC 05 Nov 1983) in the Marge series was 6 h after the time of 
recurvature. Table 4 gives the symbols used on the storm charts in the remainder of the 


study. 


TABLE 4 
STORM CHART LEGEND 


Best track positions (12 h intervals) 
Current position @ 
12 h old position es 
Future positions e 


Forecast track positions (12 h intervals) * 


Positive vorticity contours 


Negative vorticity contours 





Each chart is analyzed by first checking around the storm center location to 
determine the dominant contour pattern. In Fig. 3.3a, the first reported position of 
Typhoon Marge 1s located on the southern edge of a closed zero contour at 400 mb. 
The major pattern in the contours 1s a primary node in the zero and -2 contours, which 
has an orientation of approximately 310 degrees. Secondary nodes are found at 340 
deg in the 0 contour and 325 deg in the -2 contour. Three distinct centers are found in 
the subtropical ridge, which is well to the north of the storm at this time. Notice that 


lines from the storm center through the primary and secondary nodes appear to pass 


24 


through the breaks between the centers in the subtropical ridge. The vorticity contours 
on the northern side of the ridge appear to have little relationship to the tropical 
cyclone as they are influenced more by the positive vorticity centers associated with 


two extratropical cyclones to the north of 40°N. 


TABLES 


BEST TRACK - RELATIVE VORTICITY PATTERN 
ANGULAR CORRELATION EVALUATION CRITERIA 


Included Angle (a) Quality Assigned Points 


Excellent 6 
ood 5 

Pair 

Poor 

None 





To describe the correlation between the contour pattern and the best track, a 
comparison of the angle between the contour feature and the best track position is 
made. A point value is assigned according to the magnitude of the angle between the 
pattern axis and the best track. Table 5 gives the breakdown of points and a qualitative 
description for the ranges of the included angle (a) considered. 

In applving the evaluation scheme to Fig. 3.3a, the closed zero contour that 
the center lies on is not considered. The orientation of the primary node in the open 
zero contour is evaluated as excellent (6 points) at 42 h and the secondary node is 
judged to be fair (4 points) at 42 h. These times are estimated by drawing an arc from 
the storm position and the maximum extent of the contour pattern to the best track, 
and interpolating along the best track at the intersection point (if the time is greater 
than 72 h, a value of 72 is assigned). For example, the primary and secondary nodes on 
the -2 contour are evaluated as excellent and good respectively at 72 h. 

The vorticity contours in the 700 mb and 400:700 mb laver average fields (Fig. 
3.36 and c) have less detail near the storm, although both have a node to the 
northwest in the zero contours. The 700 mb node in the Zero contour is rated excellent 
at 58 h. The small +2 closed contour is evaluated as excellent at 12 h. The laver 


average (LAV) node in the zero contour is rated to be excellent at 48 h. Notice that a 
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Fig. 3.3a Relative vorticity (x10°s"!) at 400 mb at 12 UTC 30 Oct 1983. 
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distinct break in the ridge is defined by the -2 contour north of the 72 h best track 
position. 

The movement and development of these basic patterns is followed through 
the next eight days in 12 h intervals to illustrate the progression of these synoptic 
pattern-track relationships. The objective is to determine if the synoptic pattern 
changes precede storm track changes consistently, which would allow them to be useful 
for forecasting the track change. 

As development occurs in the next 12 h, a new +2 contour is found around 
the storm (Fig. 3.4a). Although the orientation of the zero contour nodes change little 
in 12 hours, the -2 node has moved slightly south and closer to the storm center in 
response to a strengthening of the subtropical ridge to the northwest. Notice that the 
change in the -2 node correlates with the direction change in the best track at 60-72 h. 
The strenghtening of the ridge to the northeast seems to be linked to the passage of a 
vorticity center (45°N, 170°E) moving to the east-northeast, with strengthening of the 
ridge following passage of the +12 center and weakening in proximity of the +6 
center northwest of Marge. The 700 mb field (Fig. 3.4b) has changed little, other than 
a widening of the ridge in the meridional direction and the enclosure of the storm 
center by a +2 contour. Similarly, the LAV (Fig. 3.4c) has an enclosed vorticity center. 
Although the ridge has strengthened, the break in the mdge near 25°N, 140°E has 
closed. This again suggests the storm may not continue the northerly track. 

The primary (secondary) node on the Zero contour at 400 mb becomes more 
(less) pronounced in the next 12 h (Fig. 3.5a). However, a large secondary node 
develops on the -2 contour. Although the new node on the -2 contour is larger than 
the one that has been tracked previously, it is considered as secondary due to the lack 
of past history for the feature. Furthermore, the new node does not correlate with the 
primary zero contour node or the storm track history well. At 700 mb (Fig. 3.5b), an 
expansion of the vorticity center occurs around the tropical cyclone, a closed zero 
contour develops and some weakening of the ridge occurs. A distinct change has taken 
place in the LAV chart (Fig. 3.5c) as a primary node in the zero contour has developed 
to the northwest of the storm center. This node is rated as a good correlation (5 
points) at 60 h. However, the subtropical ridge has also widened slightly and a break 
has redeveloped to the north of the storm center in response to the mid-latitude trough 
pattern. The forecaster would have to evaluate whether this break in the ridge would 
signal continued recurvature, or whether the development of a strong primary node to 


the northwest would indicate a turn to the west. 
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Fig. 3.4a Relative vorticity (x10°s74) at 400 mb at 00 UTC 31 Oct 1983. 
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Fig. 3.4b Relative vorticity (x10°s7!) at 700 mb at 00 UTC 31 Oct 1983. 
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Fig. 3.4c Relative vorticity (x10>s7!) in LAY at 00 LTC 31 Oct 1983. 
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Fig. 3.5a Relative vorticity (x10°s74y at 400 mb at 12 UTC 31 Oct 1983. 
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Fig. 3.5b Relative vorticity (x10°s7) at 700 mb at 12 UTC 31 Oct 1983. 





Fig. 3.5¢ Relative vorticity (x10°s"!) in LAV at 12 UTC 31 Oct 1983. 
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At 00 UTC O1 Nov 83, the 400 mb relative vorticity field (Fig. 3.6a) has slight 
changes including a +2 closed contour around the tropical storm center and 
strengthening of the ridge to the north of the center. Although continued growth of the 
secondary node in the -2 contour is observed, the tropical cyclone appears to be 
embedded in the primary node that in this case is associated with sustained motion 
toward the northwest. In Fig. 3.6b, the continued expansion of the vorticity center of 
the storm is evident as well as the appearance of a pronounced node in the zero 
contour to the northwest. The large vorticity gradient oriented NW to SE to the 
northwest of the storm would seem to suggest continued movement to the northwest. 
The LAV pattern (Fig. 3.6c) also shows slight changes. Although the subtropical ridge 
break has closed, a definite pinched ridge pattern remains, and the meridional 
separation of the -2 contours is reduced. The zero contour has changed little anda +2 
contour has enclosed the center. The -2 contour changes would seem to indicate that 
recurvature should be anticipated at this trme. However, the zero contour at all levels 
continues to display a major axis in the northwest-southeast direction, which would 
indicate movement to the northwest. 

At 12 UTC 1 Nov 83, a significant change occurs at 400 mb (Fig. 3.7a) as the 
zero contour now has a pronounced secondary node simular to the node in the -2 
contour. Thus, the direction of the major axis of the pattern is less clear than before. 
However, these nodes are quite far north of the tropical cyclone. Also, note that the 
major axis of the +2 closed contour around the center remains in the northwest- 
southeast orientation. Based on the persistence of the primary node and the recent 
track of the storm, no change is made in the primary node designation at this time. 
The 700 mb chart (Fig. 3.7b) shows a direct link to the higher latitudes has developed 
in the form of a narrow channel through the ridge near 22°N, 130°E. The LAV chart 
(Fig.3.7c) has a more conservative change in the contour pattern with only a gradual 
growth of the secondarv node on the zero contour. Although the ridge to the north of 
the center is quite broad, the ridge to the west-northwest appears to be narrowing in 
response to a positive vorticity center now over Honshu and moving to the east. This 
center eventually will create the path to the mid-latitudes that STY Marge will follow 
in recurvature. In hindsight, the key to the present situation is to recognize that the 
opportunity for recurvature into the weakness in the ridge will not occur. Only a feint 
to the north actually occurred, presumably because the primary vorticity channel 
oriented NW-SE was a dominant factor rather than the possible influence of the more 
distant secondary 
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Fig. 3.6a Relative vorticity (x102s7h) at 400 mb at 00 UTC 01 Nov 1983. 
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Fig. 3.66 Relative vorticity (x10>s7!) at 700 mb at 00 UTC O1 Nov 1983. 
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Fig. 3.6c Relative vorticity (x10°s7!) in LAV at 00 UTC O1 Nov 1983. 
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Fig. 3.7a Relative vorticity (x10°s7!) at 400 mb at 12 UTC 01 Nov 1983. 


O° QE ee 





Fig. 3.76 Relative vorticity (x10°s"4) at 700 mb at 12 UTC OL Nov 1983. 
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Fig. 3.7c Relative vorticity (x10°s7!) in LAV at 12 UTC O1 Nov 1983. 
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The next chart sequence (00 UTC 02 Nov 83) shows a return to a definite 
northwest track orientation in the 400 mb and LAV charts. The sharp secondary node 
is no longer present in the zero contour at 400 mb (Fig. 3.8a). The 700 mb chart (Fig. 
3.86) shows the further development of the channel pattern near 20°N, 130°E, which 
connects with the +2 vorticity center south of Honshu. Although a clear opening in 
the ridge is evident to the northwest in the LAV (Fig. 3.8c), the small positive center 
detected on the previous LAV chart has moved to the south of Honshu. Nevertheless, 
a weak cell in the subtropical ndge remains directly between the mid-latitude trough 
and the tropical cyclone. 

At 12 UTC 02 Nov 83, strong indications of recurvature are seen in all three 
charts (Fig. 3.9). In Fig. 3.9a (400 mb), the primary node on the zero contour 1s clearly 
to the west-norihwest. This primary node on the zero contour shows excellent 
correlation to the best track (6 points) at 36 h. However, a definite shift 1s seen in the 
primary node on the -2 contour in response to the continued strengthening of the 
positive vorticity center near 32°N, 142°E (good rating 5 points at 72 h). This change 
in the -2 contour nearly breaks the ridge at 24°N, 135°E. By contrast, the secondary 
node on the -2 contour shows little change at this tume. Recurvature 1s expected late in 
the period, which 1s confirmed by the good (5 points) rating at 72 h. The cyclone center 
is in a channel of positive vorticity at 700 mb (Fig. 3.9b), with the channel wall (zero 
contour) to the north nearly describing the eventual recurvature path 66 h before the 
recurvature takes place. The LAV (Fig. 3.9c) indicates the impending recurvature in 
another way. With the eastward retreat of the -2 contour that previously separated the 
positive vorticity center to the north from the tropical cyclone center, the distance 
between the zero contours separating the centers is reduced as the positive center south 
of Honshu moves to the east-northeast. Although a clear pathway to the mid-latitudes 
appeared to develop in this case, the early portion of the track was toward the west- 
northwest along the primarv node in the zero contour. 

The changes in the individual levels (Fig. 3.10a and b) at 00 UTC 03 Nov 83 
are somewhat confusing. At 400 mb (Fig. 3.10a), the ridge strengthens substantially, 
which leaves only a narrow break to the northeast of the storm center that is aligned 
with the secondary node in the zero contour. A reduction of the primary node in the 
zero contour is also evident. The channel at 700 mb (Fig. 3.10b) that had persisted for 
36 h has closed, which leaves a weakness in the ridge similar to that at 400 mb. The 


LAV chart (Fig. 3.10c) reflects the alignment of weaknesses at the upper and lower 
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Fig. 3.8b Relative vorticity (x10°s7!) at 700 mb at 00 UTC 02 Nov 1983. 
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Fig. 3.9a Relative vorticity (x10°s"!) at 400 mb at 12 UTC 02 Nov 1983. 
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Fig. 3.9b Relative vorticity (x10°s"!) at 700 mb at 12 UTC 02 Nov 1983. 
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Fig. 3.9¢ Relative vorticity (x10°s"!) in LAV at 12 UTC 02 Nov 1983. 
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bounds of the layer, and also indicates a flattening of the primary node on the 0 
contour, as well as an increase in the secondary node amplitude. 

All three charts at 12 UTC 03 Nov 83 show strong recurvature tendencies as 
the indications of 24 h previous are reconfirmed. The 400 mb and LAV charts (Fig. 
3.11 a and c) both show substantial increases in the node to the northeast of the 
cyclone center. At 400 mb, the primary node still shows excellent correlation with the 
best track (6 points) but the time factor has been reduced to about 30 h. Not only have 
the secondary nodes in the 0 and -2 contours increased in amplitude, they are aligned 
with each other along an axis that is nearly parallel with the post-recurvature track of 
STY Marge. The LAV chart (Fig. 3.11c) shows a greater flattening of the primary node 
in the zero contour (rated as excellent or 6 points at 36 h) and a marked increase in the 
amplitude of the secondary node. A line drawn from the cyclone center through the 
secondary node and the greatest southward extent of the zero contour to the north of 
the ridge also correlates well with the post-recurvature storm track. At 700 mb (Fig. 
3.116), the return of the southwest-northeast channel defined by the zero contours is 
clearly indicated. The axis of the channel correlates well with the patterns in the other 
two layers and the post-recurvature track. Consequently, the key question at this stage 
is not Whether recurvature will occur, but when? At the present speed and direction of 
the storm, it would arrive at the flattened portion of the nearest zero contour in about 
30 h. A turn toward the north should be expected at that time. 

At 00 UTC 4 Nov 83, the designation of a new primary node must be made 
on all three charts. Although the present track of the cyclone is still aligned with the 
old primary node, the positive vorticity area around the cyclone center has a distinct 
shift in the orientation of the major axis toward the northeast. Thus, the old secondary 
node becomes the new primary node. Although the 400 mb and LAV charts (Fig. 3.12a 
and c) show similar patterns with large nodes to the northeast, an additional pattern 1s 
also present. The -4 contour on the 400 mb and the -2 contour on the LAV display a 
pinched ridge pattern (Fig. 3.2d) that aligns well with the post-recurvature path. The 
700 mb chart (Fig. 3.126) includes an extended channel pattern far to the northeast 
that encompasses the “wake” of the positive vorticity center that started as a small 
center south of Korea (Fig. 3.6c) and tracked to the east-northeast to the north of STY 
Marge. The low center passed to the north of STY Marge at 12 UTC 02 Nov 83 (Fig. 
3.9c) and passed from the chart area after 12 UTC 03 Nov 83. That is, a frontal trough 
remains at 700 mb along the expected recurvature path. This feature splits the 
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Fig. 3.10a Relative vorticity (x10°s7!) at 400 mb at 00 UTC 03 Nov 1983. 
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Fig. 3.10b Relative vorticity (x10°s"!) at 700 mb at 00 UTC 03 Nov 1983, 
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Fig. 3.10c Relative vorticity (x10°s"!y in LAV at 00 UTC 03 Nov 1983. 
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Fig. 3.11a Relative vorticity (x10>s7!) at 400 mb at 12 UTC 03 Nov 1983. 
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Fig. 3.11b Relative vorticity (x10°s"4) at 700 mb at 12 UTC 03 Nov 1983. 
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Fig. 3.llc Relative vorticity (x10°s7}) in LAV at 12 UTC 03 Nov 1983. 
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subtropical ridge around 25°N and the mid-latitude ridge that extends southward from 
dO°’N, 100°E. 

By 12 UTC 04 Nov 83, an extensive node of positive vorticity has penetrated 
into the subtropical ridge during the past 24-36 h. In Fig. 3.13a (400 mb), the node axis 
continues to be nearly parallel to the best track after recurvature with a slight right 
bias. Thus, the recurvature track appears to punch directly through a rather substantial 
ridge at 400 mb. However, Marge is clearly approaching the ridge axis at 700 mb (Fig. 
3.13b). The zero contour around Marge at 700 mb and the one associated with the 
mid-latitude positive vorticity band are much closer than those at 400 mb or previously 
at 700 mb (Fig. 3.12b). This indicates that the subtropical ridge is narrowing at the 
lower levels more than at 400 mb at this time. Both lower and upper-level charts need 
to be examined in these baroclinic (verticallvy-sheared) conditions. The LAV chart (Fig. 
3.13c) has a combination of a broad node to the northwest and a narrowing of the 
subtropical ridge. The LAV node axis is not as clearly correlated with the post- 
recurvature track as is the 400 mb depiction. This poor correlation may be due to the 
layer averaging routine that does not use a mass weighted scheme, but rather a simple 
averaging of the two input level values. That is, the 700 mb level should have much 
more weight in the LAV fields because it represents a deeper layer of the atmosphere 
that 1s steering the storm. 

Over the next 12 h, the tiemeds continue in all three charts (Fig. 3.14 a, b and 
c). A continued reduction of the distance separating the zero contours of the ridge is 
seen. The next 12 h period is not available as the 12 UTC 05 Nov 83 fields are missing 
from the archived GBA data. Thus, the transition at the time of recurvature between 
OO and 12 UTC cannot be observed for Marge. 
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Fig. 3.12a Relative vorticity (x10°s7)) at 400 mb at 00 UTC 04 Nov 1983. 


58 





: | ; ie 
4 Z LSGe~ 160 | \ | 


Miges 2b Relative vorticity (x10°s7!) at 700 mb at 00 UTC 04 Nov 1983. 
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Fig. 3.12c Relative vorticity (x10°s"!) in LAV at 00 UTC 04 Nov 1983. 
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Fig. 3.13b Relative vorticity (x10>s7}) at 700 mb at 12 UTC 04 Nov 1983. 
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Fig. 3.14a Relative vorticity (x10°s7!) at 400 mb at 00 UTC 05 Nov 1983. 
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Fig, 3.14b Relative vorticity (x10°s7!) at 700 mb at 00 UTC 05 Nov 1983. 
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Fig. 3.14c Relative vorticity (x10°s"!) in LAV at 00 UTC 05 Nov 1983. 
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A major development of a positive vorticitv center at 400 mb (Fig. 3.15a) 
occurs over southern Japan at 00 UTC 06 Nov 83. Consequently, onlv a very narrow 
ridge separates the tropical cyclone and the mid-latitude circulation. At 700 mb (Fig. 
3.156), a +2 contour surrounds both the tropical and mid-latitude systems with a large 
node in the zero contour on the east side of a north-south channel. The +2 closed 
center in the node has excellent correlation with the cyclone’s final best track position. 
The LAV (Fig. 3.15c) includes two distinct centers in a positive vorticity channel with 
about 300 n mi separation of the +2 contours around the centers. A reflection of the 
700 mb node 1s present in the -2 contour to the east. Again, a more representative 
vertical weighting in the LAV would probably display the recurvature conditions more 
accurately. 

Over the next 24 h, the track of STY Marge appears to be determined by that 
of the mid-latitude system tracking to the east-northeast. At 400 mb (Fig. 3.16a and 
3.17a), the tropical center seems to be moving slightly faster than the mid-latitude 
system. In the 700 mb charts (Fig. 3.16b and 3.17b), only one vorticity center 1s 
depicted. Notice the extreme vorticity gradient between the tropical cyclone and the 
ridge to the east. The best track tropical cyclone position is to the southeast of the 
large center. The LAV charts show a blend of the two levels as expected (Fig. 3.16c 
aa! 7C). 
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Fig. 3.15a Relative vorticity (x10°s"!) at 400 mb at 00 UTC 06 Nov 1983. 
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Fig. 3.156 Relative vorticity (x10°s7!) at 700 mb at 00 UTC 06 Nov 1983. 
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Fig. 3.15c Relative vorticity (x10°s7!) in LAV at 00 UTC 06 Nov 1983. 
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Fig. 3.16a Relative vorticity (x10>s"!) at 400 mb at 12 UTC 06 Nov 1983. 
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Fig. 3.16b Relative vorticity (x10°s7) at 700 mb at 12 UTC 06 Nov 






| a . 
\ 


\ % : 
~~, fie Bh 

i y cus 

140 ? / 460 ' _- 


Fig. 3.16c Relative vorticity (x10°s7}) in LAWat?2 UTC 06 Nov 1983. 
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Fig. 3.17a Relative vorticity (x10°s7!) at 400 mb at 00 UTC 07 Nov 1983. 
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Fig. 3.176 Relative vorticity (x10°s7!) at 700 mb at 00 UTC 07 Nov 1983. 
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Fig. 3.17c Relative vorticity (x10°s7!) in LAV at 00 UTC 07 Nov 1983. 
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3. Summary of the Relative Vorticity Analysis of STY Marge’s Track 

The development of the node pattern in the relative vorticity contours and the 
change in the pattern with time appear to give reasonable indications of the three 
critical changes in the track of STY Marge. The first step in the track occurred at 12 
UTC 31 Oct 83 and was the smallest of the track changes. Each of the 400, 700 and 
LAV chart series indicated that the subtropical ridge was strong. The primary nodes in 
the vorticity patterns showing excellent correlation 18-30 h beyond the 20 deg left turn 
in the track. The strongest indicator was the 400 mb pattern, which provided the best 
indication of the conditions during the first step in the storm track. 

The weakness in the ridge that accompanied the bending of the storm track to 
the right after 12 UTC on the 3lst is evident in the 400 mb and LAY charts, but the 
pattern does not appear clearly at 700 mb. The ridge at 700 mb narrows to some 
degree but remains intact to the north of the cyclone center. The ridge break that 
develops to the northeast of Luzon 1s well to the west-northwest of the center. In this 
case, the pattern at 700 mb would indicate west-northwest movement rather than 
recurvature, about 48 h prior to the left turn at 12 UTC on 2 Nov 83. The 400 mb 
chart 24 h prior to the step in the storm track (Fig. 3.7a) could lead to a serious track 
forecast error if it was used without the other two charts. The sharp north-northeast 
node would seem to indicate that recurvature should take place at about 140°E. 
However, the 700 mb chart indicates a relatively strong ridge at 140°E with a weakness 
in the ridge to the west near 130°E. The average of the ridge weaknesses at the two 
levels would be at 135°E, which coincidentally would have been an excellent forecast as 
STY Marge eventually did recurve at about 134°E. 

The critical recurvature period appears to be well represented in the relative 
vorticity fields as the 400 mb node pattern and the 700 mb channel pattern began to 
align in the vertical after the second track step. The transition of the primary node 
pattern from one oriented toward the NW to become oriented with the previous 
secondary node started 72 h prior to the recurvature, and was nearly complete 36-48 h 
prior to the turn to the northeast. This would seem to provide a useful forecast aid if 
such patterns appeared consistently. Finally, orientation of the node and channel axes 
gave an indication of the post-recurvature track. The question is whether the patterns 
observed in the relative vorticity fields of STY Marge are evident in other tropical 
cyclones. Table 6 gives the dominant vorticity pattern, rating value and time factor for 
each 12 h period for STY Marge. 
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D. APPLICATION OF THE ANALYSIS TECHNIQUE TO OTHER STORMS 
1. Dominant Vorticity Patterns 

The relative vorticity contour pattern analysis described above is applied to 
the other 11 storms in the data set. It is not possible to describe these cases in the 
detail given above for STY Marge. Dominant and secondary contour patterns (Fig. 
3.2) are identified and the angular difference (Table 5) between the pattern and the best 
track is measured, and the time interval over which the technique would provide useful 
guidance is estimated as above. Table 7 is a summary of the relative vorticity contour 


patterns observed in the storms. 


UAE) sam| 
RELATIVE VORTICITY CONTOUR PATTERNS OBSERVED 


Storm Dominant Patterns Secondary Patterns 
Marge Node Ridge deformation, 
hanne 

Bess Node, Ridge deformation, 
Channel Resultant vector 

Gordon Node, Ridge deformation, 
Channel Storm-storm (Fave) 

Forrest Closed centers, 


Ridge deformation, 
Node 


Resultant vector 


Lex Node Channel 
Ida Node Closed centers, 
Resultant vector 
Clara Ridge deformation, Closed centers. . 
Node Storm-storm (Bill) 
Thad Ridge deformation, Closed centers 
Channel 
Vanessa Node, Storm-storm (Thad, Warren) 
Channel 
Abby Channel, Ridge deformation, 
Node Closed centers. 
Storm-storm (Ben, Carmen) 
Resultant vector 
Wayne Node, 
Channel 
Ike Channel, Node 


Closed centers 


The node pattern was found in 11 of the 12 tropical cyclones. Only Typhoon 
Thad did not show the node pattern. Thad was unusual in that it displayed a broken 
ridge from the start. The break gave an excellent indication of the 72 h position in Fig. 
3.18. Official forecast errors were large due to numerical model predictions for 
strengthening of the ridge, which did not verify (Joint Typhoon Warning Center, 1984). 
The channel pattern (Fig. 3.2a) was the next most frequently occurring pattern. 
Channels were dominant in portions of the tracks of seven storms and secondary 
patterns in two other tropical cyclones. They were generally found in the low level field 
(700 mb) or in the late stages of a storm. The mdge deformation pattern (node 
intrusion, pinch, or break) was a dominant pattern in three storms (including Thad) 
and as a secondary pattern in four storms. The closed center pattern occurred almost 
as often as the ridge deformation pattern, but was a dominant pattern only in Typhoon 
Ike. The track of Ike (Fig. 3.19) is within the channel south of the zero contour toward 
the closed +2 center before [ke turned toward the node to the west of Luzon. 

2. Cyclone-Cyclone Interaction 

The interaction between tropical cyclones is a complex problem for the 
forecaster. A determination first has to be made whether or not an interaction is 
occurring. The key questions are at what distance between centers will the interaction 
begin and when will this interaction be strong enough to change either the track or 
intensity of one or both of the cyclones. According to Sandgathe (1987), storm-storm 
interaction takes place in four of the systems in the data set. STY Vanessa interacted 
to some degree with two other tropical cyclones, but the track of Vanessa (not shown) 
does not appear to have been influenced during these interactions. The first interaction 
was a reduction of upper-level outflow above the incipient Vanessa due to deformation 
of the upper-level field by the outflow of TY Thad (Joint Typhoon Warning Center, 
1984). During the second interaction, Vanessa influenced the track of TY Warren, 
which was located in the South China Sea and had been tracking slowly westward. As 
STY Vanessa approached recurvature in the Philippine Sea, Warren turned rapidly to 
the east-northeast and increased in translation speed. This created a major forecast 
problem because of the potential danger to a group of ships leaving Subic Bay. 
Another storm-storm interaction between TY Clara and TY Bill resulted in little or no 
change in the track of Clara, but had a profound effect on Bill’s track (not shown). 
The circulation and intensity of Clara at 12 UTC 18 Nov 1984 (Fig. 3.20) were very 
large compared to that of Bill. Clara’s recurvature track was not changed by Bill, but 


Bill had a long track to the southeast under the influence of the circulation of Clara. 
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Figs 3S” Relative vorticitv (x10°s7}) in LAV at 12 UTC 18 Oct 1984 
for Typhoon Thad. Forecast tracks are dashed with stars. 
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Fig. 3.19 Relative vorticity (x10°s"!) in LAV at 12 UTC 30 Aug 1984 
for Tvphoon Ike. 
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Fig. 3.20 Streamline analysis at 925 mb at 12 UTC I8 Nov 1984 
illustrating Tvphoon Bill and Tvphoon Clara interaction 
(Joint Typhoon Warning Center, 1984). 


Typhoon Gordon interacted with TY Fave at 00 UTC 30 Aug 82 (Fig. 3.21). 
The interaction is unusual because of the nearly 14 deg of latitude (840 n mi) distance 
at which it occurs. This distance seems particularly large considering that Fave is a 
very small cyclone. The interaction may be enhanced bv the large positive vorticity 
(low) center northwest of Japan. 

In the fourth case of cyclone-cvclone interaction, STY Abby interacts with 
two other cvclones at extended distances, and apparently before the other svstems are 
even in Warning status. Usually when cvclone-cvclone interaction 1s considered, two 
mature or nearly mature systems that are approximately 600 n mi apart are involved. 
STY Abby displavs a very different pattern in the relative vorticity fields. About 36 h 
before the first warning on Tropical Storm (TS) Ben, an apparent interaction takes 
place between STY Abby (110 kt) and the tropical disturbance that would become TS 
Ben. The zero contour in the LAV field for 00 UTC 11 Aug 83 (Fig. 3.22a) surrounds 


Abbv and a +2 center to the east from which Ben will develop. The zero contour also 
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Fig. 3.21 Relative vorticity (x10°s7!) in LAV at 00 UTC 30 Aug 1982 
illustrating Typhoon Gordon and Tvphoon Faye interaction. 
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Fig. 3.22a Relative vorticity (x10°s"!) in LAV at 00 UTC I1 Aug 1983 
illustrating Super Typhoon Abby - Typhoon Ben - Typhoon Carmen interaction. 
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Fig. 3.22b Relative vorticity (x10°s"!) in LAV at 12 UTC Il Aug 1983 


illustrating Super Typhoon Abby - Typhoon Ben - Typhoon Carmen interaction. 
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Fig. 3.22c Relative vorticity (x10°s"!) in LAV at 00 UTC 12 Aug 1983 
illustrating Super Typhoon Abby - Typhoon Ben - Typhoon Carmen interaction. 
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extends toward a +2 center in the South China Sea that will spawn Tropical Storm 
Carmen 3 days hence. The pattern 1s further complicated by the strong -2 center to the 
north of Abby. In the next 12 h (Fig. 3.22b), the interaction between Abby and the 
disturbance in the South China Sea becomes stronger at a distance of approximately 15 
deg of latitude (900 n mi). The first available position for TS Ben is at 00 UTC 12 Aug 
(Fig. 3.22c) near 24.5°N, 145°E. In this field, interaction is taking place over a major 
portion of the western North Pacific basin and the South China Sea. The interaction 
may not just be multiple cyclone as large positive vorticity centers already exists that 
are adjacent (but separate from) both Ben and Carmen. The apparent interaction 1s 
taking place at a distance of 19 deg (1140 n mi) from the South China Sea to Abby 
and 16 deg (960 n mi) from Abby to the center to the northeast of Ben. Again, the 
strong ridge pattern north of STY Abby with the large negative (anticyclonic) center 
northeast of Ben further complicates the forecast situation. 
3. Forecast Performance 

Another consideration is which contours are giving the best overall forecast 
patterns and at what levels are the best forecast results found. Among the subjective 
aspects of the proposed technique are that different vorticity contours and different 
pressure levels may be involved as the storm evolves. The average rating values and 
time factors for the data set, by storm, are given in Table 8. Table 9 is a summary (by 
pressure level and by vorticity contour) of the occurrences, the average of the rating 
values and the average time factor (described above). Secondary node values are not 
given for the 700 mb level due to the small number of occurrences. Secondary nodes on 
the zero contour do not give high average ratings in this sample. The LAV field 
provides slightly better indications of track changes in this data set. This field may do 
even better if a mass weighted averaging scheme was used, as discussed previously. The 
-2 contour provides the most consistent information at all levels when it is present. 
Unfortunately, this contour is not as common as the zero contour in this data set. 

A test of the technique performance 1s how well it predicts the observed storm 
track turns. Table 10 is a compilation of the turn indication performance (hindcasts in 
this developmental sample) of the relative vorticity pattern technique. With the 
exception of the period of recurvature and northwest movement for STY Abby, the 
technique appears very successful in this sample. Typhoon Gordon also had a vorticity 
pattern with a stationary negative vorticity center to the east of Japan that is larger 


than the one that caused some problems in the forecast of STY Abby. In the case of 
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TABLE 9 
RELATIVE VORTICITY CONTOUR PERPOR IAS Grab jenn 


Level Contour # of Fests Avg Rating Avg Time 
(mb) (points) (h) 
400 2 77 4.5 (good) 35 
0(p) 169 4.4 (fair) 50 
O(s) 78 2.9 (poor) 55 
2 87 4.8 (good) 58 
700 2 42 4.7 (good) 48 
0 169 4.3 (fair) 49 
-2 62 4.5 (good) 58 
LAV 2 60 4.7 (good) 38 
0(p) 172 4.6 (good) 49 
O(s) 67 1.9 (poor) 44 
-2 69 5.1 (good) 61 


(p) primary, (s) secondary 





TABLE 10 


TURN PREDICTION PERFORMANCE 
OF THE RELATIVE VORTIG@IT PAG TERN TECRENIQOUE 


* STY Abbv - Resultant-vector pattern suggested recurvature while 
primary node indicated northwest movement. 


| 

| Turn type Observed Forecast False Indicators 

| Right 8 ? () : 
| Lett 4 4 0 

| Recurvature i 7 * 

Loop fu Q 0 





Gordon, recurvature would be forecast based on this technique. However, the pattern 


Was not as clear as in the other storms, and the time factor was reduced. A clear 
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recurvature pattern did not appear until 24 h before the event. The experience with 
this set suggests that a blocking high is a definite complicating factor in forecasting 
recurvature. 

The one right turn that is not indicated by the synoptic pattern technique is a 
30 deg change in the track of TY Ike as the system turned toward the west under the 
subtropical ridge after moving to the southwest for approximately two days. This 
occurred at low latitudes (9°N) under conditions where recurvature was not a 
consideration. 

The loops in the the track of STY Abby and of TY Lex were short-term 
features (12 h or less), and were not evident in the 12 h time step of the fields. In any 
case, forecasting loops would probably cause the Typhoon Duty Officer’s qualifications 


to be reevaluated and is not a high priority. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


This study indicates that relative vorticity charts have the potential to be a useful 
tool in the recurvature forecast decision process. Following the poor performance in 
the official forecasts of Super Typhoon Marge during October-November 1983, the 
JTWC initiated an internal review of the forecast process. This review did not reveal 
anvthing that could have been done with the existing maps and resources to improve 
the performance! This study of Marge suggests that a different presentation of the 
data available in the JTWC during the forecast development might have contributed to 
better forecasts. Thus, the first recommendation 1s that software be made available to 
JTWC to allow presentation of the Global Band Analysis fields in terms of relative 
vorticity contours. These vorticity fields allow the relative strengths and weaknesses of 
the major troughs and ridges to be analyzed readily, and incorporate dynamical 
features that are not apparent in the streamline presentation. A systematic procedure 
of locating the strong or weak points in the ridge may refine the recurvature longitude 
selection process. A mass-weighted LAV is suggested based on the experience with 
Marge. 

The developmental study had plotted best tracks on the relative vorticity contour 
charts to identify and locate the patterns described above. A true test of the 
operational importance and forecast capability of the technique should be carried out. 
The test should present an independent sample of storm sets in a manner simular to 
that found in the forecast center. The data presented should include only products and 
fields that would be available to the operational forecaster. Ideallv, a first forecast 
should be developed for each 12 h period without having the relative vorticity fields 
available, followed at a later time by a second set of forecasts developed by the same 
forecaster naving the relative vorticity fields. The skill level for each of the forecasts for 
a data set could be compared to a climatologv-persistance forecast. The use of the 
same forecaster for each set would preclude individual inherent skills from biasing the 
results. To reduce any bias that would be introduced by the forecaster remembering 
tracks, a method of presentation that mixes the sequence of storms would be required. 
If the method helps eliminate from consideration forecast aids that are not likely to be 


correct, the translation speed may be better defined by the remaining aids. Further 


IZ 


study may also reveal a pattern in the relative vorticity fields that will help determine 
the translation speed. With the proper preparation and planning, the test could be 
carried out on a quasi real-time basis utilizing the data collection and display 
capabilities in the IDEA lab. 

The storm-storm interaction analysis as currently practiced in the forecast center 
could be refined through the use of the relative vorticity fields. The pattern seen in 
STY Abby and in the TY Gordon-TY Faye cases suggest that interaction may be 
taking place at greater than expected ranges and that complex interactions patterns 
may be involved. The ability to detect these features in the analyzed fields may give 
the forecaster a basis for evaluating guidance from a numerical model that considers 
only a single vortex in the field. The early detection of vortex interactions might alert 
the forecaster to these complex interactions. Thus, a study focused on the storm-storm 
interaction cases using the relative vorticity fields is recommended. 

The previous empirical orthogonal function (EOF) representations of storm- 
synoptic data by- Wilson (1984), Schott (1985) and Weniger (1987) also might be 
applied to relative vorticity fields for eventual inclusion in the decision tree forecasting 
tool. A possible approach to stratification of the data set would be compositing of the 
data by the storm-synoptic parameters. The storm-synoptic approach could lead to 
improved forecast aid selection guidance. Another approach might be to composite 
the fields in a manner similar to the approach used by Allen (1984) in the COSMOS 
technique. Application of relative vorticity information to COSMOS tvpe field 
categories could lead to the development of an improved COSMOS technique. 

Future application and testing of this technique would be facilitated if the 
database in the IDEA lab included the best track data and the archived forecast aids. 
A study could then be undertaken to compare the relative vorticity patterns with the 
forecast aids for different categories of tropical cyclones. The objective is to select the 
“best” aid based on a set of rules developed from the past performance of the various 
aids under the conditions presented. The rules might be incorporated into an artificial 
intelligence shell on JITWC’s computers for track aid selection. Even experienced 
forecasters could benefit from the detailed analysis of the new data required for 
inclusion in the rule base. Once the basic system is developed, the “experience” of the 
system could be continually updated as new information is available. This would 
provide a system for incorporating lessons learned in a particular forecast experience, 


which is now often lost when the forecaster leaves JTWC. Another advantage in the 
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development of an AI technique is that an interactive mode could be used in the 
training of new forecasters. Rules and decision points would be displayed to the trainee 
as a scenario develops. Appropriate and immediate feedback as decisions are made 
would develop the trainee’s skills in a number of situations. 

Even though new dynamical models or other track guidance tools may be 
developed, the typhoon forecaster still must evaluate this guidance in relation to 
existing techniques. This study suggests that the relative vorticity fields would be a 
useful tool in displaying and understanding the interaction of the typhoon with 
adjacent synoptic features. It is expected that availability of these vorticity fields could 


assist the forecaster and hopefully improve track forecasts in the future. 
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APPENDIX A 
PROGRAM--GLOBAL BAND FILE IDENTIFICATION 


//SHGBID84 JOB (4697,9999), SHERMAN 2200°,CLASS=G 
[/*MAIN LINES = 40 


{im 
it 


[22 GBFID NEDN 


{7% 
Ai 


ii EXEC FORTVCLG,PARM.FORT=’‘LANGLVL(66)’ 
|/FORT.SYSIN DD * 
C 
C PGM TO IDENTIFY DATA FIELDS PRESENT IN GLOBAL BAND ANALS 
C ON FNOC TAPE 144X49 MERCATOR GRID NEDN FORMAT 
E 
C ORIGINAL AUTHOR UNKNOWN, UPDATED AND MODIFIED BY 
c B. T. SHERMAN, APRIL 1987 
€ 
C THE PROGRAM AS PRESENTED WILL GENERATE A LIST OF THE DATA 
C FIELDS PRESENT ON A NEDN TAPE OF GBA FORMAT GRIDDED DATA. 
C TWO CHANGES ARE REQUIRED TO RUN THE PROGRAM. THE JOB NAME 
C AND USER ID INFO WILL NEED TO BE CHANGED TO GIVE A UNIQUE 
C IDENTIFIER TO THE OUTPUT. THE TAPE TO BE READ MUST BE 
C NAMED IN THE UNIT 10 SPECIFICATIONS AT THE END OF THE 
C PROGRAM 

REAL*4 ID(4745),BCD(23) 

INTEGER“4 IDTG(8),NP(8),ISKP(33), ITITLE(42),NC(3) 

INTEGER*4 IN1(23), ITAU(4), BCDTAU(4) 

INTEGER IDATA(7056) 

INTEGER*2 JDATA(14112) 


€ 
DATA ENDS/‘AAAA;/ 
C 
EQUIVALENCE (IDATA(1),JDATA(1)) 
C 
NA = 14220 
NR= 7056 
NB=0 
NCNT=0 
NSK =0 
C 
CALL BUFFER(10,1,14220,ID) 
C 
C SKIP TO DESIRED DATA NSKPRC=NO. OF RECORDS SKIPPED 


9D 


C 


READ (5,1) NSKPRC 

| FORMAT (15) 
WRITE (6,2) NSKPRC 

2 FORMAT (’ NO. OF RECORDS SKIPPED = ’,I5) 
IF (NSKPRC.EQ.0) GO TO 55 
IF (NSKPRC.LE.1000) GO TO 5 
NSK = NSKPRC/1000 
DO 4 MSK=1,NSK 
CALL SKPREC(10,1000,&99,&999) 

4 CONTINUE 
NCNT = NSKPRC 
NSKPRC = NSKPRC-(NSK*1000) 

5 CALL SKPREC(10,NSKPRC,&99,&999) 


SCONE 


€ 


DO 926 L=1,4745 


926 ID(L)= ENDS 


C 


3] 


CALL GETREC(10,&99 ,&9999) 
NCNT=NCNT + 1 


CALL GBYTES(ID(1),NC,0,6,0,3) 
DO 31J = 1,3 
INI(J) = NC() 
CALL GBYTES(ID(1),ITAU,24,6,0,4) 
CALL XBCD(4,ITAU,BCDTAU) 
CALL GBYTES(ID(2),IDTG,16,6,0,8) 
I=! 
DO 32J = 411 
IN1(J) = IDTG(D 
CALL XBCD(1,IN1I(J),BCD(J)) 
l=1+1 


*« CONVERT BCD TO DECIMAL NUMBER ***** 


CALL XINT(BCD,IYMDH,IER) 
CALL GBYTES(ID(4) ,NP.0.6,0,8) 
I=1 
DO 333 = 12,19 
INI(J) = NP(1) 

l=I+1 
CALL GBYTES(ID(6),IUNIT,8,5,0,0) 
IN1(20) = IUNIT 
CALL GBYTES(ID(6),ISCLE, 13,6,0,1) 
IN1(21) = ISCLE 
CALL GBYTES(ID(7),M,24,12,0,0) 
IN1(22) = M 
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CALL GBYTES(ID(8),N,4,12,0,0) 
IN1(23) = N 
CALL GBYTES(ID(27), IDATA,8,16,0,NR) 
IF (IER.NE.0) GO TO 55 
CALL XBCD(20.IN1,BCD) 
WRITE(6,203) IYMDH,(BCD(K),K = 1,20), NCNT 
203 FORMAT(1X,,-YMDH,BCD =’,1X,110,2X,3A 1,2X,8A1,2X,9Al, 
12X,’RECORD NO. =‘,I8) 
GOTO 55 
9999 WRITE(06,120) 
120 FORMAT(IX,’ PARITY ERROR GETREC ’) 
STOP 
99 WRITE(06,121) 
121 FORMAT(IX, ’ EOF ENCOUNTERED IN GETREC’) 
STOP 
999 WRITE(06,122) 
122 FORMAT(1X,’ SKPREC FAILURE”) 
STOP 
END 
C E&ERREEEEEELEREEEEEEREEEEEREEEEELEEEEEELEEEEECLEKLEEEKE 
e 
SUBROUTINE XINT(BCD,IYMDH,IER) 
REAL BCD(23),DIGITS(10)/’0’,'1’,2’,'3’,'4’,'5’,'6,, 
wT" 89" 
INTEGER INTGRS(10)/0,1,2,3,4,5,6,7,8,9/, INTS(8) 
DO 100 1=4,11 
DO 200 J=1,10 
IF(BCD(I).NE.DIGITS(J))GO TO 200 
Il=I-4+1 
INTS(II) = INTGRS(J) 
GO TO 100 
200 CONTINUE 
WRITE(6,1000) 
1000 FORMAT(2H_,,ERROR,CHECK YOUR SUBROUTINE XINT’) 
IER=! 
RETURN 
100 CONTINUE 
1Y MDH = INTS(1)*10000000 
1 +INTS(2)*1000000 
2. +INTS(3)*100000 
3. +INTS(4)*10000 
4  +INTS(5)*1000 
5  +INTS(6)*100 
6  +INTS(7)*10 
7 +INTS(8) 
IER=0 
RETURN 


on 


EMI 
C &&&&&&&&&&&&&KKKEKEKKKELLLKEEKEREEKLERELKLKKKKEKEKLKELEEKEKEEK 


SUBROUTINE XBCD(N,INPUT,BCDOUT) 


SUBROUTINE XBCD CONVERTS N INPUT CHARACTERS FROM CDC 
EXTERNAL BCD TO IBM EBCDIC. THE INPUT CHARACTERS ARE 
STORED IN ARRAY INPUT AND THE OUTPUT (CONVERTED) 
CHARACTERS ARE STORED IN ARRAY BCDOUT. 


oo @ 2h ee 2 eee 


INTEGER HEX(47)/Z31,Z32,Z33,Z34,Z35,Z36,Z237,Z38,239, 
* 221,222,£23,L.25;L20,4.20;4.2) LLnL eee 
* Z16,Z17,Z18,Z19,Z0A,Z01,Z02,Z03,Z04,Z05,Z06,Z07,Z08, 
* Z09,230,220,210,Z11,; 233,228,226 20 Dez OZ Oa 0C |, 
"el NEMEEGN ) 
REAL CHARI) A’, By CoD eg ele 


TOK LEM IN OP OUR SCT UV" 
WXYZ 0% y 2134 '5./6/7'8, 
Lh gh ea tee tan Ne ye rae ay 
* BCDOUT(N) © 
DO 100 I=1,N 
DO 200 J=1,47 


IF(INPUT(I).EQ.HEX(J)) GO TO 150 
200 CONTINUE 

J=39 
150 BCDOUT(I) = CHAR(J) 
100 CONTINUE 


RETURN 
END 

[3 

/LKED.SYSLIB DD 

// DD 

if DD 

/ DD 

/ DD 


// DD DISP=SHR,DSN= SYS3.TAPE.LOAD,LABEL =(,,,IN) 
//GO.FTLOFOO01 DD UNIT = 3400-5, VOL = SER = GB1984, 


iI DISP=(OLD,PASS),LABEL =(1,NL,,IN), 
[| DCB =(RECFM = F,BLKSIZE = 14220,DEN = 4) 
//GO.SYSIN DD * 

0 


is 
[| 
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APPENDIX B 
PROGRAM--GET GLOBAL BAND ANALYSIS 


C THIS PROGRAM IS USED TO CREATE DATA TAPES FOR USE IN THE 

C IDEA LAB FROM NEDN FORMAT TAPES OF THE GLOBAL BAND ANAL. 
C THE PROGRAM READS THE TAPE, CONVERTS THE NEDN FORMAT TO 
C IBM EBCDIC FORMAT FOR PROCESSING AND THEN WRITES THE 

C SELECTED FIELD DATA TO TAPE IN ASCII FORMAT. THE OUTPUT 

C FROM THE GBFID PROGRAM IS USED TO SELECT THE NUMBER OF 

C RECORDS TO BE SKIPPED TO ARRIVE AT THE START OF THE FIELDS 
C THAT COMPRISE A STORM DATA SET AND IDENTIFY THE RECORDS 

C THAT MAKE UP THE DATA SET. IF RECORDS ARE MISSING FROM THE 
C GBA TAPES ENSURE THAT THE DATE/TIME INFOR FOR THE MISSING 
C FIELDS ARE NOT INCLUDED IN THE LIST OF DATA INPUT AT THE 

C END OF THE PROGRAM. SPECIFIC INFORMATION ON THE PARAMETERS 
C THAT WILL HAVE TO BE CHANGED FOR EACH DATA RUN IS SHOWN 
C AT THE END OF THE PROGRAM JCL STATEMENTS, NEAR THE LINES 

C TO BE CHANGED. REMOVE ALL COMMENT LINES THAT ARE NOT IN 
C THE PROGRAM OR SUBROUTINES BEFORE ATTEMPTING TO RUN THIS 
C PROGRAM. 

@ ORIGINAL AUTHOR UNKNOWN 

C UPDATED AND MODIFIED BY B. T. SHERMAN, MAY 1987 

C 

LiSG JOB (__,9999),’ 2200’ ,CLASS=G 

/*MAIN LINES =40,RINGCHK = NO 


j/% 


i 
ae 
// EXEC FORTVCLG,PARM.FORT =’LANGLVL(66)’ 
WEOKIT.SYSIN DD * 
C 
C PGM TO READ DATA FROM FNOC TAPE 144X%49 MERCATOR GRID NEDN 
C FORMAT 
o 
DIMENSION GLAT(49),GLON(144) 
REAL*4 ID(4745),BCD(23) 
REAL*4 DATA(7056) 
INTEGER*4 IDTG(8),NP(8), ISK P(33), ITITLE(42), NC(3) 
INTEGER“4 IN1(23),ITAU(4), BCDTAU(4), TSTTWO/’2’/ 
INTEGER IDATA(7056) 
INTEGER*2 JDATA(14112) 
CS" \ O ees 
C DIMENSIONS OF OUT MUST MATCH THE GRID SIZE OF OUTPUT 


REAL OUT(144,41),D(144,49) 


7 


REAL XL(4) 


DATA ENDS/‘AAAA’/ 
DATA XL’A EG 

C A-SFC E-700 G-400 T-250 1-200 
DATA NLVL/4; 
DATA VRU/'1'/,VRV/‘2'/ 


DATA C/.0436610743/,DD/114.5915590262/ 
NAMELIST/NGRID/ NWST,NEST,NSTH,NNTH 


EQUIVALENCE (IDATA(1),JDATA(1)) 


NA = 28220 
NR= 7056 
NB=0 
INNS = 0 
NCNT=0 
INR=0 
IOVER=0 
IER = 0 
A = 100. 
B = -100. 
CCC 
C 
C THESE ARE CURRENTLY FIXED, BUT CAN BE VARIABLE 
G 
CLAT = .| 
SEON = 1739 


~-) 
CY 


DEFINE VALUES OF LAT,LON OF ALL GRID POINTS 


DO 2000 I= 1,49 
GLON(I) = 60. + (I-1)*2.5 
2000 CONTINUE 
DO 2001 I= 50,144 
GLON(I) = ((I-49.)*2.5)-180. 
2001 CONTINUE 
DO 2002 II= 1,49 
[= 50-11 
GLAT(I) = 90.-ATAN(EXP(-ABS(II-31.)*C))*DD 
IF (11.GT.31) GLAT(I)=-GLAT(I) 
2002 CONTINUE 


eon) Orman) (ae () Ole) ©) OO) 


100 


ce REAUIN GRID SIZE 


C 
READ (5,NGRID) 
WRITE (6.NGRID) 
NI= NWST+NEST+1 
NJ=NSTH+NNTH +1 
C 
IUCOMP= 20 
IVCOMP= 20 
NRECU =0 
NRECV=0 
c 
C** READ IN NO. OF TAPES TO BE USED 
¢€ 


READ (5,78) NTAPE 
78 FORMAT (12) 
CC DO 1000 ITP=1,/7 
teNIT = 10 
CALL BUFFER(JUNIT,1,14220,1D) 
Cc 
e7* READ IN THE NO. OF TIME PERIODS TO BE INCLUDED IN THE 
S STEP OF SKIPPING RECORDS AND SKIP TO DESIRED DATA 
C 
READ (5,79) KTP 
79 FORMAT (13) 
READ (5,7) NSKPRC 
7 FORMAT (15) 
CALL SKIP (JUNIT,NSKPRC,IJEND) 
MTP =0 
IYMDH=0 
224 MTP=MITIP+1 
C 
(Geter MEAL As see fe oe ahs ahs 2h he 
C 
READ (5,225,END = 123) [LY,MTHJDJH,IREW 
225 FORMAT (12,2X,312,1X,11) 
IF (IREW.EQ.1) CALL REWIND (JUNIT) 
KYM DH = 1Y*1000000 + MTH*10000 + JD*100 + JH 
IF (LYMDH.EQ.KYMDH) GO TO 75 


IF (MTP.LE.KTP) GO TO 751 

READ (5,79) LTP 

KTP=KTP+LTP 

READ (5,7) NSKPRC 

NSKPRC = NSKPRC-2 

CALL SKIP (JUNIT,NSKPRC,IEND) 
IF (IEND.NE.0) GO TO 88 
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C 


751 CONTINUE 
927 DO 926 L= 1.4745 
926 ID(L)=ENDS 


C 


CALL GETREC(JUNIT,&88,&9999) 
GO TO 331 


88 IF (NTAPE.EQ.1) GO TO 99 


CALL REWIND(JUNIT) 
JUNIT=JUNIT+1 

CALL BUFFER (JUNIT,1,14220,ID) 
READ (5,79) LTP 

KTP=KTP+LTP 

READ (5,7) NSKPRC 

CALL SKIP (JUNIT,NSKPRC,IEND) 
GO TO 75! 


331 CONTINUE 


C 
@ 


CALL GBYTES(TD(1),NC,076,0;3} 
DO 31J = 1,3 


31 INI) = NC(J) 


Cc 


CALL GBYTES(ID(1), ITAU,24,6,0,4) 
CALL XBCD(4,ITAU,BCDTAU) 


C CHECK TO SEE IF IT ISA FCST FIELD IF SOR Gi Or i 


IF (BCDTAU(4) .EQ. TSTTWO) GO TO 751 
CALL GBYTES(ID(2),IDTG,16,6,0,8) 

I=1 

DO 32J = 4.11 

INK(J) = IDTG(I) 

CALL XBCD(I,IN1(J),BCD(J)) 

=I+1 


** CONVERT BCD TO DECIMAL NUMBER *7**2" 


CALL XINT(BCD,IYMDH,IER) 

IF (IER.NE.0) GO TO 751 

CALL GBYTES(ID(4) ,NP,0,6,0,8) 
I=] 

DO 33J = 12,19 

INWJ) = NP(I) 

=l+] 

CALL GBYTES(ID(6),IUNIT,8,5,0,0) 
IN1(20) = IUNIT 

CALL GBYTES(ID(6),ISCLE,13,6,0,1) 
IN1(21) = ISCLE 

CALL GBYTES(ID(7),M,24,12,0,0) 
IN1(22) = M 
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Le aa ss Ss 


“= a ttienien —EE—E—E—————_—— 


CALL GBYTES(ID(8),N,4,12,0,0) 
IN1(23) = N 
CALL GBYTES(ID(27),IDATA,8,16,0,NR) 
c 
C CDC HAS ONE’S COMPLEMENT, IBM HAS TWO’S COMPLEMENT FOR 
eee NEGATIVES. 
DO 301 = 1,NR 
30 IF (IDATA(I) .LT. 0) IDATA(I) = IDATA(I) + 1 
C 
ee CHECK YYMVIIDDHH TO GEF THE DATA WE NEED ****** 
CALL XBCD(20,IN1,BCD) 
IF (1YMDH-KYMDH) 751,75,224 
C AE NS HS AS TES HS AS Ss AS STG AS OE SIS AS NE OS OES AS OES BES NS NS Ee Ne fe aie Ee 20s AS STS Me ahs ole oe fs he as ake ole oes ahs als ofs ate af ate sis fs aie 
75 CONTINUE 
IF(BCD(3).NE.VRV .AND. BCD(3).NE.VRU) GOTO 751 
C 
CSET UP TO UNPACK 16 BIT DATA POINTS. THUS JDATA WHICH IS 
C *2 HAS ONE DATA POINT PER WORD. DATA IS IN THE RIGHTMOST 
C 16 BITS OF IDATA WHICH IS *4. THE SCHEME BELOW DOES NOT 
C WORK IF THE NBTISDP = 20 OR 08. 
c 
L=0 
NR2= NR*2 
{F CISCLE .LT. 32) GO TO 40 


DO = R2? 
L=L+1 
IDUM =JDATA(I) 
17 DATA(L) = FLOAT(IDUM)/FLOAT(ISL) 
GO TO 45 
40 ISL =2**(32-ISCLE) 
DO 18 1=2,NR2,2 
L=L+1 
IDUM =JDATA(I) 
18 DATA(L)= FLOAT(IDUM)*FLOAT(ISL) 
45 CONTINUE 
C 
C TEST FOR LEVEL 
E 
LVL=0 
DO 500 K=1,NLVL 
LVL=K 
500 IF(BCD(1).EQ.XL(K)) GOTO 505 
501 GOTO 751 
505 CONTINUE 
€ 
510 CONTINUE 
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NONE 


QO GO 


SEGON TINUE 
DETERMINE GRID POINTS ISJEJSJE FOR [EERE tT Ora. 


CALL FINDGR(GLAT,49,CLAT, 1, NSTH,NNTH.JS,JE) 
CALL FINDGR(GLON,144,CLON,-1,NWST,NEST,IS, IE) 
IF (IS.LE.0) IS=IS+ 144 

IF (IE.GT.144) IE=IE-144 


CREATE THE OUTPUT VECTOR Obl 
IF JS.LE.0 OR JE.GT.49, DATA AT THAT GRID POINT OF 
VECTOR OUT WILL BE CONSIDERED MISSING 


JD=JS-1 
DO 550 J=1,NJ 
IDD=IS-1 
IFLAG=0 
JD=JD+1 
IF (JD.GT.JE) GO TO 551 
IF (JD.GT.0.AND.JD.LE.49) GO TO 539 
IFLAG=1 
539 DO 540 I=1,NI 
IDD=IDD+1 
IF (IDD.GT.144) IDD=IDD-144 
IF (IDD.GT.IE) GO TO 555 
IF (IFLAG.NE.0) GO TO 538 
OUT(I,J)=D(IDD,JD) 
IF (OLT(I,J).LT.B.OR.OUT(I,J).GT.A) IER = 1 
GO TO 540 
538 OLT(I,J)= 99999. 
IOVER = 1 
540 CONTINUE 
550 CONTINUE 
GO TO 570 
551 WRITE (6,552) JS,JE,JD 
552 FORMAT (‘ ERROR IN IDENTIFYING LAST LATITUDE LINE‘/’ 
1JS=',13,2X,'JE=’,13,2X,'JD = ’,13) 
STOP 
555 WRITE (6,556) IS,IE,IDD 
556 FORMAT (‘ ERROR IN IDENTIFYING LAST LONGITUDE LINE’/’ 
LIS = ,13,2X,’IE= ’,13,2X,'IDD = ‘,I3) 
STOP 
570 CONTINUE 


IF(BCD(3).NE.VRU)GO TO 800 
NRECU=NRECU +1 
WRITE(6,10) NRECU,IYMDH,(BCD(K),K = 1,20) 
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l ,CLAT,CLON IS, IE,JS,JE 
l ,D(IS,JS), D(IE,JE),OUT(1,1),OUT(NI,NJ) 
2 JER,IOVER 
10 FORMAT (1X,’NRECU =’,14,110,1X,3A1,1X.8A1 
eorex CLAT= F5.1, CLON= ,F6.1,2X. 
21S =',13, IE=",13,’ JS=",13,’ JE=‘,13,4F7.1,13,11) 
CC WRITE (IUCOMP) NRECU\(BCD(K),K = 1,3), 1YMDH,OUT 
WRITE(IUCOMP,610)NRECU,(BCD(K),K = 1,3), 1Y MDH 
610 FORMAT(I4,1X,(3(A 1,1 X)),110) 
WRITE(IUCOMP,620)OUT 
620 FORMAT(10F8.2) 
IOVER = 0 
IER = 0 
GO TO 751 
800 IF(BCD(3).NE.VRV)GO TO 751 
NRECV=NRECV +1 
WRITE(6,11) NRECV,IYMDH,(BCD(K),K = 1,20) 


1 .CLAT,CLON,IS,IE,JS,JE 
1 ,DUIS,JS),DUIE,JE),OUT(1,1),OUT(NINJ) 
2 JER,IOVER 


11 FORMAT (1X,’NRECV = ‘,14,110,1X,3A1,1X,8A 1 
1,1X,9A1,2X,’CLAT=’,F5.1,, CLON = ’,F6.1,2X, 
21S = "13, IE=',13,’ JS=',13,’ JE=',13,4F7.1,13,11) 
CC WRITE (IVCOMP) NRECV.(BCD(K),K = 1,3), |YMDH,OUT 
WRITE(IVCOMP,630)NRECV,(BCD(K),K = 1,3), 1Y MDH 
630 FORMAT(I4,1X,(3(A1,1 X)),110) 
WRITE(IVCOMP,640)OUT 
640 FORMAT(10F8.2) 
IOVER = 0 
IER = 0 
GO TO 751 
CC STOP 
C 
775 WRITE (6,776) LYMDH,IX2,NCNT 
716 FORMAT(‘ERROR IN RECORD ARRANGEMENT/7/‘IYMDH = ’,110,2X, 
l'IX2=’,110, RECORD NO.=’‘,15) 
STOP 
777 WRITE(6,150) lYMDH 
150 FORMAT(1X,,PGM ENDED IYMDH =’,I10) 
STOP 
9999 WRITE(06,120) 
120 FORMAT(1X,’ PARITY ERROR GETREC ’) 
STOP 
99 WRITE(06,121) 
121 FORMAT(1X, ‘ EOF ENCOUNTERED IN GETREC ’) 
STOP 
999 WRITE(06,122) 
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122 FORMAT(IX,’ SKPREC FAILURE ’) 
1235ShOP 
END 
(CR See a oe eee eras OO cee tet as SUBROUTTINE SKIP ME 2 Ie 2 ae Oe ae MR Oe Bh oh of He He 2% ais a ie He Me HK 
e 
SUBROUTINE SKIP GUNIT,NSKPRC,IEND) 
C 
C THIS SUBROUTINE SKIPS NSKPRC RECORDS FROM JUNIT 
S 
IEND=0 
G 


WRITE(06,2) NSKPRC 
2 FORMAT (0°*,NO. OF RECORDS SKIPPED= ’,I5) 
IF (NSKPRC.EQ.0) GO TO 6 
IF (NSKPRC.LE.1000) GO TO 5 
NSK = NSKPRC/ 1000 
DO 4 MSK=1,NSK 
CALL SKPREC(JUNIT,1000,&99,&999) 
4 CONTINUE 
NSKPRC = NSKPRC-(NSK*1000) 
5 CALL SKPREC(JUNIT,NSKPRC,&99,&999) 
6 CONTINUE 
GO lO £25 
99 IEND=1 
RETURN 
999 WRITE(06,122) 
122 FORMAT(IX, SKPREC FAILURE ’) 
STOP 
123 RERURN 
END 
C FR ESRF E SS Ets SUBROUTINE FINDGR aft fe ste ofe ofe oie te ae fs che fe Ee ote aft fe ofe fe ake ote fe of ots fe oe fs oe ok 
G 
SUBROUTINE FINDGR (A,N,P,IFLAG,NGI,NG2,MS,ME) 


THIS SUBROUTINE DETERMINES THE GRID LINE CLOSEST TO 
THE CYCLONE CENTER 
AND RETURNS THE BOUNDARY GRID LINE NOS. 


DEFINITION OF PARAMETERS 

A= ARRAY OF N VALUES OF LAT/LON 

P= POSITION (LAT/LON) OF CYCLONE CENTER 

NGI=NO. OF GRID POINTS ON WES@/ SOL Teesepe OF T 

NG2= NO. OF GRID POINTS ON EAST/NORTH SIDE OF P 

MS= GRID LINE NO. ON WEST/SOUTH SIDE OF P 

ME=GRID LINE NO. ON EAST/NORTH SIDE OF P 

IN THE CASE THAT P IS EXACTLY IN BETWEEN IWO GREIBSBIENES: 
IFLAG DETERMINES WHICH GRID LINE WILL BE GONSIDERED Wie 


te) (Gann) GiiG@iaGont.) () (a6) ©) 
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CEOSES | 

IPEAG = 1 GRID LINE Wilt HIGHER VALUE THAN P WILL BE 
Grose N 

IFLAG=-1 GRID LINE WITH LOWER VALUE THAN P WILL BE 
GolO sem 


DIMENSION A(N) 


ee) (2) ee 


IF (P.GE.-177.5) GO TO 25 

Ll = 49 
L2 = 50 
RETURN 

25 DO 30 1=1,N 
Il=1+1 
IF (11.GT.N) If=1 
IF (P.LT.A(1).OR.P.GT.A(II)) GO TO 30 
D1 =P-A(I) 
D2=A(II)-P 
IF (D1-D2) 31,32,33 

31 MS=I-NGI1 
ME=I+NG2 
GO TO 50 

32 IF (IFLAG.GT.0) GO TO 33 
GO TO 31 

33 MS=II-NGI 
ME=II+NG2 
GO TO 50 

30 CONTINUE 
WRITE (6,39) P 

39 FORMAT(‘FAILURE TO FIND LAT/LON CLOSEST TO CYCLONE 

1ICENTER’/’ LAT/LON =’,F8.2) 

STOP 

50 RETURN 
END 

C RGSS BKK SKK SKE SKE EESEREKKS 
@ 

SUBROUTINE XINT(BCD,IYMDH,IER) 


c 

REAL BCD(23), DIGITS(10)/0’,'1’,'2’,'3','4','5',6,, 

789") 

INTEGER INTGRS(10)/0,1,2,3,4,5,6,7,8,9/,INTS(8) 
C 

IER=0 
C 

DO 100 1=4,11 

DO 200 J= 1,10 


IF(BCD(I).NE.DIGITS(J))GO TO 200 
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Il=I-4+1 
INTS(II) = INTGRS(J) 
GO TO 100 
200 CONTINUE 
WRITE(6,1000) 
1000 FORMAT(2H_ ,, ERROR,CHECK YOUR SUBROUTINE XINT’) 
IER=1 
RETURN 
100 CONTINUE 
IY MDH = INTS(1)*10000000 
1 +INTS(2)*1000000 
2. +INTS(3)*100000 
3. +INTS(4)*10000 
4 +INTS(5)*1000 
5 +INTS(6)*100 
6  +INTS(7)*10 
7 +INTS(8) 
RETURN 
END 


REEKEKKEKRKEEEEEEERKEEKEEKEKKEEEEKEEKEEKKEKEEEEKEKEEERLEKEEERKE 


SUBROUTINE XBCD(N,INPUT,BCDOUT) 


C 

c 

c 

C SUBROUTINE XBCD CONVERTS N INPUT CHARACTERS FROM CDC 
C EXTERNAL BCD TO IBM EBCDIC. THE INPUT CHARACTERS ARE 

C STORED IN ARRAY INPUT AND THE OUTPUT (CONVERTED) 

C CHARACTERS ARE STORED IN ARRAY BCDOUT. 

G 


INTEGER HEX(47)/Z31,Z32,Z33,234,Z35,Z36,237,238,239, 
* 721 ,Z23,Z24,Z25,Z26,Z27,Z28,Z29,Z12,Z13,Z14,Z15,Z16, 
* Z17,Z18,Z19,Z0A,Z01,Z02,Z03,Z04,Z05,Z06,Z07,Z08,Z09, 
* Z30,Z20,Z10,Z11,Z3B,Z2B,Z2C,Z1B,Z1C,Z0B,Z0C/,INPUT(N) 
REAL CHAR(47)/‘A’BY CD’, BE’, FG Ek, 1)" 
KL OMEN OUP OURS TCU VW’, 
XY TY TED SAS 6 T ERG 
Qe af ot one ye vag rae oly"! BEDOUT(N) 
DO 100 [=1,N 
DO 200 J=1,47 
IF(INPUT(I).EQ.HEX(J)) GO TO 150 
200 CONTINUE 
J=39 
150 BCDOUT(I)= CHAR(J) 
100 CONTINUE 
RETURN 
END 


% 


CCCC ALL OF THESE LIBRARIES ARE REQUIRED FOR THE PROGR A 
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CCC WORK, THE PROGRAM ALSO CALLS A USER INACCESSIBLE ROUTINE 
een PPE DaGe TREC WHICH [IS A PART OF THE TAPE UTILITIES 

SCCee bn wer weviOVE THESE FOER EINES TO RUN* 

//LKED.SYSLIB DD 

{/ DD 

fi DD 

ii DD 

i DD 

// DD DISP=SHR,DSN=SYS3.TAPE.LOAD,LABEL =(,,,IN) 

SECC LARSEN IW OreOn FROSPAPEMENTS ARE THE READ AND 
CC WRITE FORMATS FOR THE TAPE UNITS. 

eC Tre riRST IS THE READLFORMAT 

CC FOLLOWING ‘VOL=SER=' GIVE THE TAPE NAME(ON THE OUTSIDE) 
Seerilbl IN DATA FROM THE TAPE SCAN AT THE APPROPRIATE SPOT 

eC LABEL= 

CC BLKSIZE= 

eC DEN= 

wees ReVOVeE PHeESse J COMMENT LINES TO RUN ****ee"eRe% 
/i'GO.FTIOFOOL DD UNIT= 3400-5, VOL= SER=GB1984, 

I DISP=(OLD,KEEP),LABEL=(1,NL,,IN), 

{| DCB =(RECFM = F,BLKSIZE = 14220, DEN = 4) 

CC THE NEXT CONTROL STATEMENT IS THE WRITE TO TAPE STATEMENT 
CC THAT PUTS THE DATA ON THE TAPE IN ASCII FORMAT. 

CC FIRST CHANGE THE EXTERIOR LABEL LINE ‘VOL=SER= XAAXXX 
meee. ! YOU WILL NEED TO ADJUST THE BLKSIZE OF THE OUTPUT 

Se ONE FIELD SHOULD FILL ONE BUOCK FOR EPFPICENT MANAGEMENT 
CC ON THE VAX, & BECAUSE THIS METHOD USES LESS TAPE/RECORD 


CC RECSIZE = (HEADER) + (8*GRIDSIZE) 
GC HEADER = 1*LRECL 
CC GRIDSIZE = NPTS(N/S)*NPTS(E/W) 


CC BLKSIZE MUST BE AN INTEGER MULTIPLE OF LRECL. 
CC IF RECSIZE IS NOT AN INTEGER MULTIPLE OF LRECL USE THE 
CC NEXT INTEGER TO GET BLKSIZE 
EC EX: IF RECSIZE = 63.4 
CC BLKSIZE = (64 * 80) = 5120 
CC ****REMOVE THESE 15 COMMENT LINES TO RUN***TCAD#*###* 
//GO.FT20F001 DD UNIT = 3400-5, VOL = SER = DSGBS84, 
//  DISP=(OLD,KEEP),LABEL =(1,NL.,OUT), 
// DCB=(OPTCD=Q,RECFM = FB,LRECL = 80,BLKSIZE = 28320,DEN = 4) 
eee 
//GO.SYSIN DD * 
&NGRID NWST=48, NEST=95, NSTH= 10, NNTH=30&END 
| THIS {Ss THE NUMBER OF TAPES TO BE READ 
3040 THIS THE NUMBER OF RECORDS TO SKIP 
84 69120 YEAR,MONTH,HOUR, REWIND CODE:0= NO, |= YES 
84 61000 
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APPENDIX C 
PROGRAM--GLOBAL BAND TAPE TO DISK 


program GBTTOD 


This program reads Global Band Analysis data from tape 

and writes it to a file on disk for further processing. 
The program uses QIO system services to read the data. 

The subroutine (mtread) which actually calls the QIO is 
located in diskS2:[metdept.tape]mt/lib. 

(GBTTOD means Global Band, Tape TO Disk) 


Brett Sherman & Jim Cowie 7/87 (NPS) 


KKK IKK IKK KKK HK KKH KK KKK KKK KK KKK KK KKK KKKKKKEKKKKKKKK KK KKK KKK KKKKK 


NOMA FNANANAANANAANRANO 


character data*28320 
integer*2 data_len, real len 
integer*4 istat, status, channel, sysSassign, sysSdeassign 
parameter(data_len=28320) 
c begin 
c assign the QIO channel 


istat = sysSassign( ‘tape’, channel,,) 
write(*,*)’channel #’,channel 
if(.not. istat)then 
write(*,*)’ could""nt assign QIO channel number’ 
write(*,*)’ istat=’,istat 
goto 999 
endif 


c open the disk output data file 
open (unit=1,file=’gbanal.dat’,status=’nevw’ ) 
c using QIO’s, read the record from tape and write to disk 


do 20 j=1,4 
100 call mtread(channel,%ref(data),data_ len,status,real_ len) 

write (*,*) ‘real len=’,real len 

if (.not.status) then 
call libSsignal(%val(status) ) 
go to 999 

end if 

do 10 i=1,28320, 80 

write(1,200)data(i:i+/79) 


110 


10 continue 


20 continue 
c go back and read the next header and grid 
goto 100 


c deassign the QIO channel and exit the program 
999 istat = sysSdassgn( channel ) 
200 format(a) 

stop 

end 


IT] 


APPENDIX D 
PROGRAM--FILE CONVERTER 


program FILE CONVERTER 
cc This program converts the FNOC Global Band field data located ina 
cc disk file created by GBTTOD.FOR, to a gempak grid file. the program 
cc uses gempak subroutines to create the navigation data, the grid file, 
cc & write the data to the new file. each file contains a complete 
cc storm data set. 
eG 
cc Initial data relating to the grid size and lat/lon will have to be 
cc adjusted to values corresponging to the grid being converted. Ngrid 
cc will have to be changed to reflect the number of grids in the data 
cc set. The file names for read file, header file and the gem file will 
cc have to be changed for each different storm data set. 
aie 
cc In order for the grid conversion to run properly the LAT.DAT file 
cc must be accessible by the program. The program converts the 2.5 
cc degree longitude square GBA grid to a 2.5 deg latitude by 2.5 deg 
cc longitude grid for use in the GEMPAK routines. 
aC 
cc LCDR B.T. SHERMAN & JAMES COWIE JUNE 1987 
cc 
CC-------------- -- - - = - - = - - - = - $$ - $e = 
cc 
ee 
character file*80, rfile*80, setname*80, parm*4, lvl*1, hdrfil*80 
integer level(2) 
character time(2)*13, yymmdd*6, hh*2, mm*2, lat£il*80 
logical replace 
integer fn, type, comp, grdno 
dimension grid(80,40), rnvbl1k(256) ,ighdr(2), tgrid( 7056) ,anbl1k( 128) 
real intlat,igrid(144,49),mgrid(144,40),lat(49) 


Cc 
c define the file to be read & the name of the file created 
Cc 

data rfile /’THAD 1984.DAT’/ 

data file /’THAD 1984.GEM’ / 

data hdrfil /’THAD 1984.HDR’/ 

data latfil /’LAT.DAT’ / 
 & : 

open (unit=4,file=latfil,status=’old’ ) 

open (unit=5,file=rfile,status=’old’ ) 

open (unit#6,file=hdrfil,status=’ new’ ) 
Cc 


replace = .false. 


c read the latitude data 
Cc 
jl = 1 
do 4i=1,/7 
read(4,5) (lat(j), j = j1 , jl + 7) 
5 format(8(£9.3)) 
jl = j1 + 8 
4 continue 
c do6Mm = 1,49 
ec write(6,*) lat(m) 
e 6 continue 
e 
Cc 
c set initial parameters 
€ 
c ngrid is the number of grids in the data file to be converted, kx & 


cc ky are dimensions of the grid, the lat/log limit designators are set 

cc up as follows ex: rlatb = the real value of the latitude at the 

cc bottom of the grid, mm = the minutes for the time designator, ivcord 

cc = 1 designates pressure as the vertical coordinate, navsz,ihdrsz, 

cc ianlsz are set by the subroutines in the gempak, level(2) = -1 is the 
cc indicator that you are using only one level for the grid 


~ 


ngrid = 104 

kx = 80 

igx = kx 

ky = 40 

igy = ky 

npts = kx * ky 
rlatb = -40.000 
rlatt = 57.500 
rlonr = 257.500 
rlonl = 60.000 
navsz = 256 
ianlsz = 128 
ihdrsz = 2 
iveord = 1 


data mm /’00’/ 
a 
c call the subroutine to create the nav data for the grid 
a 

call gr mllnv (kx,ky,rlatb,rlonl,rlatt,rlonr,rnvblk,iret) 
e 
c if the return code indicates a problem print it and stop 
c 

if (iret.ne.0) then 

write(*,*) ‘gr mllnv error - iret=’,iret 

go to 999 

end if 


Q 


call the subroutine to create the new data file 


call gd_cret(file,navsz,rnvblk,ianlsz,anlblk, file, ihdrsz,iflno,iret) 
if (iret.ne.0O) then 
write(*,*) ‘gd cret error - iret=’,iret 


go to 999 

end if 
2 

grdno = Q 
c 

do 300 1 = 1,ngrid 

read(5,10) f£n,lvl, type, comp, yymmdd,hh 

10 forfaat(1X,13,1x,al,1x,1i1,1x,11,3x,a6,a2) 

Cc 


WRITE(*,10)FN, LVL, TYPE, COMP, YYMMDD, HH 


grdno = grdno + 1 


write(6,15)fn,lvl, type, comp, yymmdd,hh,grdno 


ys. format(1x,i3,1x, 41, 1x%,i11,1x,1143x,aGea2piae ace) 


Q 


time(1) 
time(2) 


yynmdd//’ 4 Mii 


if (lvl.eq.’A’) then 
level(1)=1000 


else if (lvl.eq.’E’) then 


level(1)=/700 


else if (lvl.eq.’G’) then 


level(1)=400 


else if (lvl.eq.’T’) then 


level(1)=250 


end if 


level(2)=-1 


if (comp.eq.1) then 


parm =’UWND’ 


else if (comp.eq.2) then 


16 


20 
400 


parm =’VWND’ 
end if 


reaaq (5,16) 

format( 80x) 

il=l 

do 400 m=1,/706 

read (5,20) (tgrid(I), 1 =i11,11+9) 
format (10(f£8.2)) 

11=11+10 

continue 
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Cc 


Cc 


Cc 


Cc 
Cc 


Q 


Cc 


Cc 
C 


c 


write(6,*)’tgrid (1) =',tgrid(1) 


kount = 0 
do 100 5 = 1,49 
do 200 k = 1,144 
kount = kount + 1 
igrid (k,50-j) = tgrid (kount) 
200 continue 
100 continue 
interpolate the GBA grid to NOGAPS format 
do 510 i=1,144 
mgrid(i,17)=igrid(i,19) 
intlat=2.5 
intj=18 
North of the Equator 
do 520 j=20,48 
if(lat(j).lt.intlat.and.lat(j+1).gt.intlat)then 
mgrid(i,intj)=igrid(i,j)+(igrid(i, j+1)-igrid(i,j))* 
* (intlat-lat(j))/(lat(j+1)-lat(j)) 
int} = intj +1 
intlat = intlat + 2.5 
endif 
320 continue 
South of the Equator 
int} = 16 
intlat = -2.5 
do 530 j = 18,2,-1 
if (lat(j).gt.intlat.and.lat(j-1).lt.intlat)then 
mgrid(i,intj)=igrid(i,j)+(igrid(i, j-1)-igrid(i,j))* 
* (intlat-lat(j))/(lat(j-1)-lat(j)) 
int} = intj - l 
intlat = intlat - 2.5 
endif 
530 - continue 
510 continue 
make the grid to go to the GEM file 
do 540 j = 1,40 
do 550 k = 1,80 
grid(k,j) = mgrid(k,j) 
550 continue 
540 continue 
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write(*,*)grdno, time, level, parm 


call gd writ (iflno,grid,igx,igy,ighdr, time, level,ivcord, parm, 
nbsiz,replace,iret) 
if (iret.ne.Q) then 
write(*,*) ‘gd writ error - iret=’,iret 


go to 999 
end if 
300 continue 
999 stop 
end 


This latitude information is required by the grid conversion 
subroutine in the program FILE CONVERT (remove these three lines prior 
to processing data) 


-40. 
-24. 
-04. 
14. 
2a 
48. 
Be, 


956 
226 
990 
821 
986 
069 
745 


042 
-21.928 
-02.497 
1) e222 
334050 
49.711 
0.0 


-39 


~37.076 
-19.593 
00.000 
19.593 
37.076 
31.300 
0.0 


-35.056 
-17.222 
02.497 
21.928 
39.042 
526835 
0.0 
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-32.986 
-14.821 
04.990 
24.226 
40.956 
54.318 
0.0 


-30.866 
-12.393 
07.473 
26.484 
42.816 
Joe 70 
0.0 


698 
943 
943 
698 


-28. 
-09. 
GP . 
28. 
44. 
7 oleae 
0.0 


621 


-26.484 
-07.473 
12.4623 
30.866 
46.372 
58.462 
0.0 
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